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A  method  of  treatinq  polyatomic  systems,  finite  or 
extended,  is  presented  *hich  fully  exploits  their  "atcm-in- 
mclecule"  nature.  Within  an  independent-particle  model  a 
partitiorinq  technique  is  applied  to  the  projection  of  the 
full  polyatomic  space  into  many  atotric  subspaces.  The 
sutspaces  are  then  each  ccupled  to  one  another  approximately 
thrcuqh  second-order  in  overlap  in  a  piecewise  self- 
consistent  fashion.  The  inherent  localized  picture  allows 
for  an  approximate  fern;  of  the  interatomic  interactions 
without  resortinq  to  neqlecf  of  any  differential  overlap  or 
use  cf  any  etrpirical  parameters.  Molecular  orbitals  and 
enerqies  may  be  obtained  from  an  approximate  electron 
propagator  which   is  tased  on   a  ncdel  Hamiltonian   built  up 


frcir  atctric  one-electron  Hamiltcrians  in  diagonal  form. 
Symmetric  orthonormalization  of  these  orbitals  gives  a 
density  aatrix  which  can  be  used  as  a  guide  for  charge 
redistribution  within  the  system.  A  generalization  cf  the 
Hartree-Fcck  apcrox imaticn  based  or  a  statistical  ensemble 
is  employed  which  permits  the  use  cf  fractional  occupation 
numbers  in  the  atomic  configurations.  Application  cf  the 
methcd  here  includes  the  nitric  oxide  molecule  and  its 
chemisor r tion  on  a  nickel  (100)  surface. 


CHAPTER  T 
INTRODUCTION 


In  the  quantuit  mechanical  treatment  of  polyatomic 
systems,  whether  they  be  finite  cr  extended,  one  feature 
consistently  emerqes,  and  that  is  they  appear  to  be  built  up 
from  "atcms".  The  reasonable  success  of  the  cellular 
(Slater,  1934),  valence-bond  (VE)  (Gerratt,  1974),  atora-in- 
inclecule  (MP!)  (Hoffitt,  1951),  and  similar  building-blcck 
(Adams,  1971;  Gilbert,  1972)  approaches  is  indebted  tc  this 
feature.  However,  as  cue  progresses  from  atoms  to  molecules 
serious  complications  are  encountered  in  the  numerical  and 
analytic  techniques  which  have  been  so  successfully  applied 
tc  atoms.  These  complications  are  largely  due  to  the  less 
of  spherical  symmetry  and  one-center  expansions.  for 
numerical  techniques,  the  threat  of  mult idimensicral 
inteqratiens  and  quadratures  has  led  tc  cellular 
approximations  tc  the  potential  (Slater,  1934)  and  local 
approximations  *c  the  exchange  interactions  (Slater,  1971). 
Cn  the  ether  hand,  analytic  methods  are  faced  with  the 
tremendous   number  cf   conplicated  nulticenter   t fco-electron 
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integrals,  and  one  is  usually  forced  tc  sacrifice  both 
quality  and  quantity  in  the  selection  of  basis  sets  (Dunning 
and  Hay,  1977).  The  exclusion  cf  cere  electrons  by  means  of 
effective  potentials  (Kahn  et  al.,1S76)  has  become  a  popular 
leans  of  reducinq  *:he  number  of  integrals.  Pl  whole  gairut  of 
senriempirical  methods  have  also  been  devised  and  constantly 
revised  tc  help  cverccre  these  difficulties. 


One  wenders,  though,  whether  or  net  the  ability  tc  even 
describe  the  "atoms"  within  a  system  is  lest  in  all  these 
sacrifices.  This  brings  us  tc  the  theme  of  the  present 
work.  Vie  wish  to  maintain  as  much  as  possible  a  conplete 
and  accurate  description  of  the  "atoms"  in  a  system  at  the 
controlled  expense  of  an  approximate,  but  sufficiently 
riqorous,  ncrempirical  treatment  of  the  interatomic 
interactions.  Some  corollaries  to  this  theme  will  be  tc 
keep,  irtact,  the  molecular  Hamiltorian  while  partitioning 
the  discrete  eiqenspace  (the  finite  space  spanned  by  the 
discrete  eigenf unctions)  of  an  effective  cne-electron 
operator;  tc  enplcy  analytic  techniques  enabling  proper 
treatment  of  the  exchange  interactions;  to  use  suitable 
basis  sets,  such  as  Slater-type  orbitals  (STO's)  of  at  least 
"double-zeta"  (17)  quality  with  polarization  functions  (DZP) 
when  needed;  to  correctly  compute  and  retain  all  one-center 
inteqrals  and  all  two-center  one-electron  integrals;  to 
reduce  all  multicenter  two-electron   integrals  tc  two-center 
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one-electron  integrals  ty  exploiting  the  lccali2ed  picture, 
but  without  resorting  to  empiricism  or  neglect  of  any 
differential  overlap;  and  to  establish  atcmic  valency  as  a 
tetter-de finec  and  Hcrkarle  ccncep*. 


In  the  nest  chapter  the  effective  one-electron 
Hamiltonian  operator  and  the  partitioning  of  its  eigerspace 
will  le  developed.  Chapter  ITI  kill  deal  more  explicitly 
with  the  pseudopotential  which  is  derived  frotr  the 
partitioning  in  chapter  II.  The  approximate  forms  of  the 
interatomic  Cculcmbic  and  exchange  interactions  are  given  in 
chapters  IV  ard  V.  Chapter  VI  is  devoted  to  the  piecewise 
construction  of  a  model  molecular  one-electron  Hamiltonian 
which  when  used  in  the  moment  expansion  of  the  electron 
propagator  gives  rise  tc  a  model  Feck  matrix  as  its  first 
mcirent.  Solution  of  the  corresponding  eigenvalue  problem  is 
discussed  together  with  an  appropriate  population  analysis 
for  redistributing  electronic  charge.  In  chapter  VI  a  irere 
qeneral  form  of  the  Hartree-Fock  approximation  is  described 
which  allows  for  the  definition  of  a  self-consistent 
potential  frcm  a  density  corresponding  to  a  statistical 
ensemble  kith  specified  occupation  numbers.  The  problem  of 
spin-orbital  occupation  assiqnments  is  the  topic  of  chapter 
VIII.  In  chapter  IX  *he  method  is  applied  to  nitric  oxide 
and  its  c hemisor ption  or  a  nickel  (100)  surface.  Separate 
treatments  for  the  nitric  oxide  molecule,   the  nickel  uetal, 


■4- 


and  the  nickel  ( 1  GO)  surface  are  included  to  emphasi2e  the 
applicability  cf  the  method  to  various  systems.  The  final 
chapter  is  a  qensral  discussion  of  some  pros  and  cons  cf  the 
method  arii  scire  ideas  for  future  *cik  to  alleviate  any 
shortcomings  the  method  might  have. 


CHAPTEF  II 

THE  EFFECTIVE  CNE-FtFCTBCN  CPFRSTCR  ANt 
THE  PARTITIONING  OF  ITS  ETGENSFACE 


In  order  tc  illustrate  the  partitioning  of  a  polyatomic 
system  into  many  atomic  subspaces  it  is  sufficient,  at  least 
for  the  ocment,  to  consider  a  diatomic  molecule,  A-B.  In 
much  of  vhat  is  about,  tc  fcllcw,  specific  attention  will  be 
focused  cr  the  case  of  a  diatomic  trclecule.  This  is  dene 
strictly  for  clarity,  and  generalization  to  polyatomic 
systems  shall  be  made  at  various  times.  The  method  is  by  no 
means  restric+ed  tc  diatomic  molecules,  even  though  its 
generalization  miqht  not  appear  obvious  at  times.  In 
connection  with  this,  we  will  mostly  be  concentrating  on  a 
single  atom  (A)  at  any  one  time,  while  the  remaining  atoms 
(E)  will  he  considered  as  its  environment.  We  shall  usually 
denote  this  by  the  use  of  superscripts  A  or  B,  respectively; 
but  when  no  specific  reference  is  made,  atom  A  is  tc  be 
assumed.  For  cases  other  than  diatomics,  primes  will  also 
be  used.  See  Figure  1  for  other  notation  conventions  tc  be 
adhered  tc  throughout  this  and  remaining  chapters. 
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Fiaure  1.    Electronic  and  nuclear  coordinates  for  the 
diatomic  irolecule  A-B. 


The  averaqc  value  of  the  ncnrelati vist ic  many-election 
flaniltonian,  with  respect  to  a  proper  choice  of  density 
operator  (see  chapters  VT  and  VII) ,  for  the  diatomic 
molecul<=  can  te  written  in  second  quantization  as 

(A)  A       +  (A)         AA       AA  +    + 

<H>    =5         b    <a   a    >+    1/2   £         (J      -K      )<a    a    a    a    > 
i  iii  ij  ijij         ijji 

(B)  B       ♦  (E)        BE       BB  +    + 

t  ^         h    <t    b    >  +    1/2   5  (J      _K       )<t    b    b    b    > 

i  iii  ij  ij      ij        i   j    j    i 

(A)      (B)       AE       AB         ♦    + 
«•    1/2  ^        £         (J      -K      )  <a    b    b    a    > 
i        i  ij      ij        ijji 

(E)      (A)        EA       BA  «■    + 

♦     1/2    5        £  (J       "K      )<fc    a    a    t    >  (IT-1) 

j        i  jiji         jiij 

where  the  a's  and  fc's  are  field  operators  corresponding  to 
spin-orti  tals  assigned  tc  atoms  A  and  B,  respectively,  and 
the  summations  have  been  restricted  accordingly.  let  us 
assume   that   the   spin-orbitals   assigned   to   atom   A   are 
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expanded      in    e       tasis   centered  cr    A      and    the      spin-or fcitals 

assiqned    to   atom    D    are    likewise  centered    on    B.       Such    a    total 

enerqy      functional      leads      to  the      effective      one-electron 

operator 

2  A    -  1  A  A 

F      =    -1/2    V    -    Z    r       +J-K 
1  1  1A  1  1 


B    -1  B  B 

Z    r      +    J    -    K 
1E  1  1 


(TI-2) 


where      the   Coulcirb      and    exchange      operators   associated      with 
center    A    are    defined    try 


A  (A)     A       A       -1       A 

J    *(D     =   £         q    <<♦     Ir       |<J>    >    4>{1) 

1        i     i   i   1 2   i 


and 


A  (A)  A   A   -1      A 

K  <MD  =5    q  <<t  |r   !0>  4>  (1) 

1       i    i   i   12     i 


(II-3) 


(11-4) 


where  the  Dirac  brackets  here  denote  integration  over  the 
coordinates  of  electron  2  and  summation  over  its  spin.  The 
q's  are  the  spin-orbital  occupation  numbers.  Analogous 
definitions  for  the  Coulomb  and  exchanqe  operators  en  atom  E 
are  assuncd.  The  s pi r-crbi ta  Is  are  no*  restricted  to  be 
eiqensol  u  tiens  *c  a  pse udoeigen value  problem, 


F  4(1)  *    E  4(1), 
1 


(II-5) 


althcuqh  they  must  be  solved  for  self-consistent ly  since  the 
operator  depends  on  them  through  the  Coulomb  and  exchange 
operators. 
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At  this  point  we  have  dore  nothing  but  invoke  the 
Hartree-Fcck  approximation  with  the  only  twist  being  the 
localized  restriction  on  the  spin- crbi tals.  The  canonical 
Hartree-Fcck  scheme  would  require  that  the  spin-cr titals  be 
symmetry-adapted ,  that  is,  invariant  (except  possihly  a 
chanqe  in  siqn)  under  the  operations  of  *he  irolecular 
symmetry  point  qroup.  This,  of  ccurse,  is  an  unnecessary, 
tut  convenient,  restriction.  Only  the  total  state  need  be 
invariant,  en-1  such  invariance  can  te  achieved  with 
equivalent  ortitals  (Hurley  et  ai,  1953) ,  which  transform 
into  each  ether  under  the  operations  of  the  group.  Such 
ortitals  are  irore  qeneral  and  can  usually  be  well-localized. 
We  shall  use  such  localized  equivalent  orbitals  in  cur 
treatment  and  restrict,  instead,  our  spin-orbitals  to  be 
centered  en  a  single  atomic  site.  So  each  "molecular"  spin- 
ortital  will  te  associated  unambiguously  to  a  specific  atom. 


We  wish  now  to   simulate  an  atom-in- molecule  partitioning 

of  the  pseudoeigenvalue   problem  in  (11-5).    To   effect  the 

partitioning,   let  us   lock  at  the  diatomic   molecule  A-B  at 

very  larqe  interatomic  separation  where   the  fifth  and  sixth 

terms  in  (TI-1)   become  negligitle  and  the  eigen functions  of 

F  in  (II-2)  approach  the  atcmic  canonical  Hartree-Fock  spin- 

crtitals  associated  with  the  two  eigenvalue  problems 

AA  2A-1A         A         A 

F  4>  |"C1  =  1-1/2  'v  -  Z  r   +  J  [C]  -  K  [0]  }  <t>  [0]    (II-6) 
i  1A  i 

and 


E    E  2B-1E  E  B 

F    <J>    r  C  1    =     (-1/2    S    -    Z    r       +    J    [C]    -    K   [ft]    }     $    [0]  (II-7) 

i  1B  i 

where  th€  f  0 1  denotes  this  sepa rated-atcm  case  (dropping  the 

reference  to  electron  1  on  th°  operators  and  the  functions). 

It   shoulc  fce   enphasized  that   this  limiting   feature  is   a 

consequence  of  the   localized  form  which  we   have  chosen  and 

that  the   canonical  molecular  crbitals  would   have  separated 

rather   unpredictably   (particularly    where   high   synmetry 

exists).    Now   as  he  allow  the   two  atoms  to   approach  each 

ether  and  interact,  let  us  focus  our  attention  on  a  specific 

atom,  say  A.   Suppose  that  instead  of  solving  (II-2)  fcr  all 

of  its  eicensclu tiers,  we  seek  only  these  solutions  which  we 

associate  with  atom   A,   while  the  solutions   of  (II-7)   for 

atom  B    are  assumed   to  suffice,   for   the  moment,    as  the 

solutions  of  (IT-2)   associated  with   atom  E.    That  is,   we 

wish  to  solve  the  pseudceiqenva lue  problem 

A  2    A  -  1    A        A 

F  *  T  1]  =  1-1/2  v  -  Z  r   ♦  J  [1  ]  -  K  [  1] 

i  1A 

E  -1    E         B         A 
-  Z  r   ♦  J  [C]  -  K[0]  J  ♦  [1]    (TI-8J 
1B  i 

for  all  the  functions  centered  en  A,  while  holding  the 
functions  on  B  fixed  as  thouqh  they  were  already  self- 
consistert  eicen functions  of  F  in  (II-2) .  The  [1]  here 
denotes  the  first  iterate  of  this  higher-level  self- 
consistent  process.  In  order  fcr  the  functions  en  tctb 
centers  A  and  B  to  re  si  irultanecus  eiqe  nf  unctions  of  F,  they 
must  fce  ncnin teract inq  throuqh  F,  that  is, 


-10- 


A        E 
i       1 


(IT-9) 


must  te  satisfied.    It  is  sufficient,  though,   that  they  be 
crthcqcral, 


<*  r  111*  f  Cl>  =  C,    ¥  i,j, 
i     i 


(11-10) 


if  either  is  in  fact  an  eiqenf uncticn.  But  before  we  treat 
this  restriction  on  the  eigenvalue  problem,  let  us  follow 
through  uith  cur  preposition  a  little  further. 


We   otviously   could   have   just    as   easily   chosen   tc 

concentrate   en  atctr   E   and  ended   up   with  the   eigenvalue 

problem 

B  2  B    -1  B  B 

F|f1l=     1-1/2    v    -    2    r       ♦    J    [1  ]    -    K   [  1  ] 
1  IB 

A    - 1  A  A  B 

-    Z    r      ♦   J    [C]   -    K   [0]   }    ♦   p]       (H-11) 

subject  tc  the  constraints 

A      B 
<♦  r0]|«>  [  11>  =  C,  f    i,i-  (11-12) 

i      1 

This   would   cive   us   a  full   set   of   functions,    <>[1],s, 

corresponding  tc  our  first  level  of  iteration.   It  is  fairly 

easy  to  see  hew  cne  could   continue  this  reasoning  to  obtain 

a  set   of  4f  2  V  s  from   the  $P"|*s,    and  so  on,    until  seme 

degree  of  convercence  is  attained,  that  is,  until 
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A  A 

♦  fnl  — >  *  r  n-1], 
i         i 

E  B 

4  r  n  1  -->  4  r  n-1  ],  (IT-13) 

1         1 

and  we  satisfy  the  crthcqcna lit y   conditions 

A      B 
<4  rnl|4  f  r.l>  =  r,    v  i,j.  (II-1U) 

i     1 


Getting   tack   now    to   the   problem   of    handling   the 

additional  constraints  en  the   pseudoeigenvalue  profclen,   we 

intend  tc   treat  this  problem   with  the   projection  operator 

technique  (LOwdin,   1961).    The   preceeding  development  and 

much  of  what  is  about  tc  follow  could  just  as  well  have  been 

treated  vith  partitioninq  technique  (LBwdin,   1964),   but  we 

shall  dispense   with  this  equivalent   approach  for   the  time 

teinc.    let  us  start  by  considerirq  a  general  unconstrained 

function  vhose  projection  is  the  desired  function, 

A    _A      (E)   3   A    B 
*  =  4   -5   <4  14  >  4 

i     i    1      j   i   j 

_A 
=  (1-P  )  4  (IT-15) 

B    i 

where 

♦     (E)  I         E 
P   =  F   =  2    14  ><*  I  (IT-16) 

E     E    j       j    j 

is  the  projection   operator  which  projects  cut   the  "B-part" 
of  the  functions  on  A.    Condition   (11-14)   in  terms  of  the 
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unconstrained    functions       (dropping    row      the    [n]-notaticn      at 

self -consistency)     becomes   automatically    satisfied, 

A  B  _A  +  E 

<<J>     |F|*    >    =    <$     |  (1-F    )     F\i)    > 

i  i  i  B  1 

_A  r  _A       +         B 

=    <$     |FU    >    -    <♦    I  P    F 1  4    > 

i  1  i      B        i 

E_AE  (E)E_AE  BB 

=    5    <<J>    |  $    >    -   ^         S    <  <f:    |  <|>    ><<*     |c»    > 

i      j  k  kikkj 

=    0  (H-17) 

where    we      have    made   use      of   the      fac+:    that    the       functions   on 

center    E    are    eiaenf unctions    of    F    and    thus   orthogonal    amongst 

themselves.         The    eigenvalue    prcblett.      associated    with    atom    A 

in    terms    cf    these    unconstrained    functions    would    te 

_A  A  _A 

F(1-F    )    4>       =    E     (1-E    )     ♦    ,  (11-18) 

E         i  i  E         i 

and  the  secular  problem  would  then  become 

_A        ♦  _A 

|<<»  M1-P  )  (F-5)(1-P  )I0  >|  =  0.  (11-19) 

i      B  B    i 

Thus,     we    have   transferred    the   restriction    on   the 

eiqenfunctions   to   a   modification  cf   the   operator.     We 

therefore   seek  instead   the   unconstrained   solutions  to   a 

£I:°JJ§cted   Hanijtonian.     The    analogous   secular   problem 

associated  with  atcti  E  would  be 

_B        +  _B 

|<<|  |  (1-P  )  (F-5)(1-P  )|«  >|  =  C  (11-20) 

i       A  A     j 
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with  analogous  definitions  fcr  the  unconstrained  functions 
on  atom  E  and  fcr  the  projection  operator  associated  with 
the  "A-part"  cf  A-B. 


If  we  rewrite  equations  (11-19)  and  (IT-20)  with  the 
unprojected  operator  F  separated  out,  we  get  (Weeks  et  al. , 
1969) 


i     A     i 


and 

_b      _e 

|<<J>  |  F  +V  -6  |*  >!  =  G 

i     D     i 

where  the  pseudcpotentials  are  qiver  by 

V      =    -  (P   F-SP    )  (2-f    )     -   [F,P    ] 
ABES  B 


and 


V      =    -  <P    F-SP    )  (2-P    )    -    [F,P     ]. 
B  A  A  A  B 


(IT-21) 


(11-22) 


(11-23) 


(11-24) 


It  should  be  noted  that  the  second  term  in  (11-23)  and 
(11-24)  will  vanish  if  the  operator  F  commutes  with  the 
prciecticr   operators,    that    is, 


tf,p  i  =  r*,p  i  =  o. 

A  B 


(11-25) 


Furthermore,    if   the   "projection"    operators   are   truly 
idetnpoteiit,  then  these  pseudcpotentials  reduce  to 


V   =  -(F-S)P 

A  E 


and 


V   =  -<F-6)P  . 

E  A 


(11-26) 


(11-27) 


-m 


As  the  problem  has  beer  formulated  sc  far,  both  of  these 
requirements  are  fulfilled,  but  these  terms  shall  be 
retained  for  reasons  which  will  be  clarified  later. 


At  this  point  the  diatomic  problem  has  been  effectively 
reduced  tc  twc  coupled  atonic  problems  without  any  serious 
approximations  outside  the  independent-particle  model  or 
limitations  imposed  on  the  basis,  but,  we  still  have  not 
really  siuplified  the  problem  with  respect  to  size  or  effort 
either.  we  have  merely  reformulated  the  problem  tc  prctrpt 
new  insiqht.  The  fact  that  v.e  have  net  yet  gained  anything 
is  evident  when  we  realize  that  the  Coulomb,  exchange  and 
proiecticr  operators  are  all  defined  in  terms  of  the 
constrained  functions  and  that  the  unconstrained  functions 
are  neither  eiqe rf urctiens  cf  F  rcr  crthonorraal.  Although 
the  unconstrained  eigenf unctions  cf  our  modified  operators 


_  A  A_  A 

(F+V  )  *  =  S  4> 

A    i  i  i 

_B  E_B 

(F+V  )  <»  =  S  i 

b   i  i  i 


(11-28) 


(11-29) 


share  a  common  set  of  eigenvalues  with  the  constrained 
eicenfunctiens  cf  F  (Weeks  et  al.  ,  1969),  they  are  not 
simultaneous  eigenf unctions  cf  the  same  operator.  Since  the 
eiqenf unctions  en  different  centers  correspond  to  different 
operators,  they  need  net  be  orthogonal.  In  fact,  because  of 
the   enercy  dependence   cf   the   pseudopotentials,   all   the 
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eiqenf unctions  should,  in  principle,  correspond  to  different 
operators.  This  enerqy  dependence  of  the  pse udopotentials 
will  be  taken  up  in  the  next  chapter. 

In  principle,  one  could  solve  directly  for  the 
unconstrained  ei  qer.  f  urct  ions  of  the  [codified  operator  and 
then  reqenerate  the  constrained  eiqenf unct ions  cf  F,  from 
which  one  could  obtain  new  Coulomb,  exchanqe  and  projection 
operators  (and  a  total  energy) .  However,  this  reintroduces 
multicenter  additions  to  our  otherwise  single-center 
unconstrained  functions,  and  thus  multicenter  integrals, 
which  is  precisely  what  we  have  set  out  to  avoid.  This 
brines  us  to  cur  first  aajcr  ap  pro  sanation. 


In  the  beqinninq,  we  assumed  that  each  spin-orbital  could 
be  unambiquous] y  assigned  to  a  particular  atom  and,  mere 
importantly,  expanded  in  a  single-center  basis.  Cur 
unconstrained  eiqenfunctions  have  this  property  imposed  on 
them  directly  by  limiting  the  basis  in  the  solution  cf  the 
projected  eigenvalue  problem.  However,  as  we  have  just 
noted  above,  this  indirectly  forces  the  constrained 
eiqenfunctions  cf  the  unprojected  eigenvalue  problem  to 
deviate  from  this  distinction.  Before  making  any  mere 
assumptions  about  the  individual  spin-orbi tals,  we  must 
recall  that  the  only  thinq  which  has  any  direct  consequence 
is  the  fcrm  of  the  total  density  itself, 
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(A)  h*    A 
$    0   0, 

i     i   i 


(11-30) 


which  occurs   in  the   Cculcmb  and   exchange  operators,  the 

projection  operators,   and  the  total  energy  expression.  In 

terms  of  the   unconstrained  functions  the  total   density  nay 
te  written 


(A)       ♦  S*      _A 
^    (1-P  )  $   (1-p  ) $  , 
i        B    i      B   i 


(IT-31) 


which  upon  expansion  of  the  projection  operators  gives 
(A)_A*_A      a)      (B)   A   B    B*_A 

5      $    o    -  5      5      <i  i*  >  o    4 

i  iiii  ijji 

(A)     (B)       B    _A      _A*    E 

-5        5        <0    |<*    >    4      0 

i         1  jiii 

_  (A)      (B)     _A       B  E    _A  B*    B 

+   ^         J         <♦    U    ><H*    >   <l      *    •  (11-32) 

i         j,k         i      i         ki         j      k 

If    we    new    invck€    the      flulliken    approximation    (1949)       for    the 

differential    overlap. 


t>    t      =    <i    I*    >(♦*    +    **)  /2, 
i   i  i      i        i   i        i   i 

then    (11-31)     becomes 


(11-33) 


-17- 


(A)_A*_A              (A)      (B)  _A       B    2  _A*_A          B*    E 

£          ♦       <t>       -    5         £  <<*     !4>><<t)       4*<*       4> 

i           iiij  ij  iijj 

(A)      (B)  _A       E    2  B*    E 

i       i  i     i  i     i 

(A)  (B)     _A  E    2       _A*_A 

i  i             i  j             i      i 

(A)     A*  A 

i           i  i 


where 

A 


(E)     _A       E    2    -1/2_A 
(>    '    =    <1    -  5         <*!<!>>)  i 

i  i  i      i  i 


and 


A         A 
<$    » | $    »>    =    1 . 

i         i 


(11-34) 


(11-35) 


(11-36) 


Thus,  the  Kulliken  approximation  leads  directly  to  a  fern  cf 
the  tctal  density  ir  terms  of  the  single-center 
unconstrained  functions  which  have  been  renor realized  tc 
unity.  So  with  this  as  an  incentive,  we  shall  replace  the 
constrained  functions,  wherever  they  occur,  by  their 
renoriralized  unconstrained  counter  j  arts ,  that  is, 


A     A 
•   =  $  '. 
i     i 


(H-37) 


We  will  drop  the  "prime"  notation  and  keep  in  mind  that  this 
substitution  shall  always  be  in  effect.  The  full  impact  of 
retaininq  the  localized  single-center  expansions  will  not  be 
realized  until  approximate  fcr&s  for  the  interatomic 
interactions  are  introduced   in  charters  IV  and  V. 
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Eefore    lookinq      in    ircre      detail    at      the    pseudcpotentials , 

let    us   ccrsider    vihat    additional      problems    we   can    expect    when 

the      sani€        reasoning      is        extended      to         molecules      teycnd 

diatomics.         In      brief,       the    secular    problem,         analogous    to 

(11-19)  ,         that      would    be      associated      with      at.otr    A      in      the 

triatoraic   molecule    A-B-E*    is 

_A  +  +  _& 

|«»    |   (1-P    )     (1-P       )      (F-S)  (1-F       )  (1-P    )   I*    >l     =0  (TI-38) 

i  B  B'  B'  B  1 

where  no*,    in  qer.eral,   the   projection  operators   are  not 
orthogonal  (or  disjoint),  that  is, 


(E)   (E«)   B   E«     E»       B 
P  P   4  O)  =  5   £     <4>  It   ><*   I4»  *  (  1) 
E  B'       i    k      j   k    k      i 


where 


E   E' 
i   k 


(11-39) 


(11-40) 


We  would  qet  analogous  expressions  for  atoms  B  and  B1.  The 
corresponding  pseudcpotentials  would  thus  be  mere 
complicated  than  those  encountered  in  the  diatomic  case, 
unless  scire  simplifying  assumptions  can  be  made.  An  attetrpt 
will  te  iade  in  the  next  chapter  to  retain  a  pseudopotential 
which  is  accurate  to  at  least  second-order  in  interatomic 
overlap.  This  next  chapter  is  devoted  to  a  more  detailed 
lock  at  the  pseudopotentials,  with  particular  attention  paid 
to  the  erergy  dependence. 


CHAPTER    III 
THE    ENEEGY-CEPFNDENT    FSEUDOPOTENTIA t 


In    the   previous    chapter,       we      have    essentially   defined    an 

effective   atom-in-mclecule      one-electron    Ha  iriltorian.  Its 

associated    eicenvalue    problem    has    the    form    us   that 

+ 
Q    FQ    t    =    S    <t>  (III-1) 

where   F  is   the  polyatomic   analog  of   (TI-2)   and   Q  is   a 

product  ct    prciecticn  operators, 


C  =  (1-P   )  d-P   )  (1-F   )  ...  . 

B       B»      B" 


(III-2) 


which   are  no^   in   general  disjoint   (and  so   Q   is  net   in 
general  a  projection  operator  itself).   If  v*e  expand  Q, 


C=1-P   -  E   -  .  .  .  +  P   P   +  F   P   +..., 

B     B«  B   E»    E   B" 


(III-3) 


and  retain   only  terms   up  to   second-order  in  differential 

overlap  jsecend-crder  ir.   P)  ,   and  then  invoke  the  Mulliken 

approximation  (1S49)    for  the   terns  involving  interatomic 

differential  cverlap, 

B     B1  B    B'     B     B      B  *    B  ' 

I*   ><*   I  =  V2<<J)   14   >{!♦   ><<t   l+!<*   ><<*   I),  (III-4) 
i     i  i    i     i    i     j     j 

we  get  the  gereral  form  for  Q 
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C  =  1  -  P 

where 

(E)        (B'J   E    B»  2    B     B 
F  =  ^    d-5    <0   !♦>)!♦  ><4   I- 

i        j#i     i    j      i    i 


(ITI-5) 


(IIT-6) 


The  first  suit  in  (III-6)  is  ever  all  normalized 
(noncrthcqonal)  unconstrained  eigenf unctions  associated  with 
all  the  atoms  (R)  in  the  molecule  except  the  one  with  which 
C  is  associated  (A),  and  the  second  inner  sum  is  over  all 
atoms  except  A  and  the  particular  atom  B  being  summed  over 
in  the  outer  sum.  The  operator  P  thus  "projects"  out  the 
space  complementary  to  the  atom  being  considered  (A),  which 
is  just  a  generalization  of  the  projection  operators  we 
encountered  for  the  diatomic  case  in  the  previous  chapter, 
that  is. 


♦      (B)  B   E    E 

p  =  p    =  S     y  u  ><*  i 

i     i   i    i 


with 


E         (B»)   E    R'  2 

y    =  1  -  5   <o  io  > 

i      i*i    i   i 


(III-7) 


(III-8) 


beinq  an  overlap  correction  to  the  "projection"  manifold. 
In  the  limit  of  zero-overlap,  F  is  a  true  projection 
operator  %ith  the  idempctency  relation  fulfilled. 


The  secular  problem  which  we  then  wish  to  solve  is  of  the 
same  forir  as  (TI-19),  namely 
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|<<*  I  (1-P)   (F-S)  (1-F)  \i>    >|  =  C 
i  i 

or  in  teens  of  a  pseudo pctentia 1 

|<$    I  F  +V-S  J<>    >|    =    0 

i  i 


(III-9) 


(111-10) 


where,    analoqcus    to    (11-23), 

V    =    -  (EF-r, t)  (2-P)    -    fF,P].  (TII-11) 

In  order  to  simplify  the  form  of  V,  as  in  (TI-26),  we  need 
to  examire  FP  and  PP  in  seine  detail.  let  us  first  consider 
the   former, 


(B)    B      B        B 

ff  =  f  5       y    i*  ><<>   I . 

i  i      i         i 


(111-12) 


He  must  recall  that  the   eigenf urctiens  in  (111-12)   are  not 
eiqenf unctiens  cf  F  but  rather 


B    B     E  B 
(F  +  V  )  <»   =  S  4 
i     i  i 


(IIT-13) 


with  the   corresponding  pseudopcten  ti  al  associated   with  the 
environment  cf  atom  B.   Rearranging  (111-13)  gives 


B    BE     BE 
F$  =  £<t>   -  V  <(> 
i    i  i       i 


(111-14) 


which    through    first-order    in    perturbation    theory    fceccmes 


B         _B    E 
F   4      =   M 
i  i    i 

where 

_E  E  BEE 

&      =    S      -    <<J>     |V     |<D    >. 
i  i  i  i 

Usinq  this  now  in  (111-12)  giv^s 


(111-15) 


(ITI-16) 
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(3)_E    REE 
EF    =   5         &    X     !♦    ><♦    I # 
i  i    i      i         i 


(111-17) 


and    since    F    is    a    sel£-ad jcint    operator,     (TTI-15)     also    grants 
FF    =    FE  (111-18) 

causinq    the   ccnmutator    in    (111-11)     to    vanish,    leaving    us    with 
V    =    -  JF-&)  (2P-PP).  (TTT-19) 

locking    new    at    FF    we    have    that 


<D    )      (E  •)      E       E  »       D  E  E'  3' 

pf  =  5       5       y    y    i*    ><<*    i <»    ><4    i . 

i  1  i      1        i  i        i  j 


(IIT-20) 
If    we    again    irvoke    Mulliken's   approximation    (ITT-U)     then 

(111-21) 


|B)_E    E       E  E 

pf  =  5       Y  Y    \i   ><i    \ 

i  i    i      i         i 


where 


_E  E  (E')     E*        B  B»     2 

y=y+5        y<cji|(*>. 

i  i        i#i         i        i        i 


(111-22) 
into    (111-19)     we   obtain 

(iir-23) 


(E)  _D  _E       E       E  E 

v  =  ^       (6-b  )  (2-y  )y  i*  ><♦  I 

i  i  i      i      i         i 


which   in    the      limit    of    an   orthonormal      "projection"    niarifcld 
reduces    tc   the    Phil  li  ps-Kleinira  r    (1959)     pseudopotential 


<E)  E  E  B 

V    =   ^  (S-6    )I4    ><*     I. 

i  i         i         i 


(ITI-24) 


Cne  should  ncte  the  explicit  dependence  of  V  cr  the 
energy  eigenvalue  S.  This  energy  dependence  is  the  crly 
reinainina   cctr  pi  ica  ticn  in   writing   dewn   the  final   matrix 
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equations  for  solvinq  the  secular  problem  in  (IIT-10). 
Eefore  qettinq  into  the  enerqy  dependence,  though,  let  us 
introduce  a  basis  and  proceed  tc  set  up  these  matrix 
eq  cations . 

For  the  sake  of  clarity  and  simplicity,  we  consider  an 
crthonornal  basis  centered  on  the  atom  (A)  *ith  which  we  are 
currently  concerned,  such  that 


4)   -  $  u  C 
i   k= 1  k  ki 


(ITI-25) 


cr   in   matrix   notation 

4    =    u    C  (IIT-26) 

where    M    is   the    number    cf    functions      in    the    basis     (and    so    the 
traxitrum    number    of      eiqenf unctions    we    can    obtain).  We    then 

consider    a   composite    basis    for      all    the    ether  eigenf unct ions 
centered    en    th^ir   respective    atomic   centers,    that    is, 


E  |B)     E    B 

*  =   ^         u    C 

i  k           k    ki 

cr  in  matrix  notation 

E  E  B 

#  =  u  c  . 

If  we  now  define  the  matrices 


(111-27) 


(ITI-28) 


S    =  <u  | u  >, 
kl      k   1 


_  ^ 


kl 


(B)    _B   BE   B* 

(2-y  )>  c  c 
i   i  ki  li 


(111-29) 


(111-30) 


and 
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(B)_B         _B       B    B       B* 

b      =  5       s   ( 2-y  )  y  c    c 

kl         i  i  iikili 


(111-31) 


then 


V        =    <u    |  V|  u    > 
kl  k  1 

t  t 

=    Z     (SES    )       -     (SBS    ) 

kl  kl 


(TIT-32) 


Thus,      aside    from      the   energy   dependence    of      V,       the    secular 
problem     (111-10)  reduces    tc    finding      the    unitary      matrix   C 

such    that 


C     (F+V)C    =    8 

where 

f         -    <u     |F|u    > 
kl  k  1 

and    fi    is    now    a    diagonal    iratrix. 


(ITI-33) 


(TIT-34) 


As  already  stated,  the  only  remaining  complication  is  the 
explicit  enerqy  dependence  of  the  pseudopotential  (TII-32). 
As  lonq  as  this  enerqy  dependence  is  there,  we  can  not  seek 
simultaneous  eiqenf unctions  of  the  same  operator,  since  each 
eiqenf unction  would  correspond  tc  a  different  operator 
dependent  upcn  its  cwn  eigenvalue.  Such  eigenf unctions 
would  net  even  be  orthogonal  tc  each  o*her.  Operators  of 
this  nature,  which  correspond  tc  a  one-dimens icna 1  Hilbert 
space,  are  alsc  inherertly  rcn-Hermitear  in  any  itatrix 
representation  (cf   order  hiqher  than   one).    We   shall  now 
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derive  ar  approximate  form  for  the  pseudopoten tial  which  is 
explicitly  enerqy  independent  and  Hermitean.  First,  let  us 
rewrite  (3H--3)  in  terns  of  a  modified  Fock  matrix,  that 
is, 


C  FC  =  S 

where 


(111-35) 


F  =  F+8P-E, 


P  =  SEE  , 


an  d 


(IIT-36) 
(111-37) 


P  =  SBS  . 
Then  we  ncte  that  from  (111-35) 


(111-38) 


FC  =  C6.  (111-39) 

No*  puttinq  (111-36)  back  into  (111-35)  we  have 

♦     ♦_     +_ 
C  FC  +  C  PCS-C  PC  =  8  (111-40) 

which    upon    substitution    of    (TI1-39)    becomes 

+  ♦_  + 

C    FC  +  C    PFC-C    PC    =    S  (III-U1) 

or 


C     (F  +  EF-B)  C    =    S. 

So  ccmparinq  this  with  (111-35)  we   have  that 


(111-42) 


F    =    F+FF-B 

which    upcr    rearranqinq    qives 

-  1 
F    =     (  1-P)        (F-B)   . 


(TlI-!43) 


(III-4U) 
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If  we  ncvi  expand  the  inverse,  we  generate  the  series 

F  =  F-F+PF-PH+PPF-.. .  (111-45) 

from  which,  alcnq  with  (111-36)  and  (111-43),  we  gather  that 

SP  =  FF  =  EF-ii+EPF-...  (111-16) 

which  is  similar  tc  a  perturbation  expansion.  Since  the 
pseudopotential  matrices  (111-3  7,111-38)  are  already  seccrd- 
crder  in  overlap,  this  series  should  be  rapidly  convergent. 
In  fact,  we  shall  take  cnly  the  first  two  terms  of  (111-46) 
in  our  approximation.  The  non-Herniiticity  is  also  apparent 
in  this  expansicn  form,  and  so  we  choose  the  approximate 
Hermitean  form  for  (ITI-46) 


SP  =  (FF-P5  +  FP-HP)  /2. 
Our  modified  Fock  matrix  then  beccnes 


(III-U7) 


F  =  F-E+(FE-BP+PF-PR) /2. 


(TTI-M8) 


We  have  at  this  point  arrived  at  a  scheme  for  treating  a 
specific  atom  within  a  polyatomic  system.  This  schene  can 
he  used  fcr  each  atom  in  the  system,  which  are  then  coupled 
in  a  self-consistent  manner  as  described  in  the  previous 
chapter.  Ml  aspects  of  the  Hartree-Fock  scheme  (or  any 
qenerali2ation  of  this  independent-particle  model)  have  been 
preserved,  with  the  environment  of  each  atom  being  reduced 
to  the  inclusion  of  additional  potential  terms.  Our  largest 
sacrifice  was  irace  in  return  for  the  one-center  expansion  of 
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the  eiqenfunctic rs;  hcwever,  we  have  not  yet  used  this  to 
its  full  advantaqe.  Althouqh  we  have  eliminated  all  three- 
and  four-center  twc-electron  integrals  and  many  two-center 
inteqrals,  we  have  only  replaced  their  by  just  as  many  cre- 
and  two-center  inteqrals.  This  reduction  alone  may  very 
well  have  beer  worth  the  sacrifice,  but  the  goal  here  is  for 
a  more  substantial  reduction  in  work.  The  localized  picture 
which  cones  ouf  of  ^his  prc-jected  Hamiltcnian  approach  can 
be  used  to  its  fullest  advantage  in  approximating  the 
interatomic  Cculcmtic  and  exchange  intpractions  which  now 
have  their  closest  resemblance  tc  an  electrostatic  ircdel. 
This  will  be  the  topic  cf  the  next  twc  chapters. 


CHAPTER  IV 

ft  PCDFL  flPERCXIMATION  TC  TEE  *FFECTIVF 
CHE-FLECTRCN  INTERATOMIC  CCtlLCMB  PCTENTIAL 


If  one  can  su tst antia lly  reduce  the  rurcber  of  multicenter 
two-election  integrals  encountered  in  conventional  analytic 
jt  initjc  methcds,  the  ccirputatiora  1  savings  would  be 
overwhelffing.  In  the  last  chapter,  the  multicenter 
inteqrals  were  all  reduced  to  at.  most  two-center  integrals 
ty  restricting  the  eiqe  r.f  unctions  to  one-center  expansions. 
The  actual  nomter  of  inteqrals,  however,  has  not  been 
chanqed.  We  shall  in  this  chapter  exploit  the  localized 
picture  and  sutstart ia 1  ly  reduce  the  number  of  two-center 
Cculcwb  integrals.  This  shall  be  done  ty  replacing  the 
effective  one-electron  operator,  for  the  two-electrcn 
Coulomb  interaction,  ty  a  ircdel  apt roxima ticn. 

The  integrals  which  we  wish  tc  approximate  are  of  the 
form 


<<t  |J  |<>  > 

i      i 


(IV-1) 
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which  is  the  Coulombic  interaction  of  the  i-th  spin-orbital 
with  all  the  otter  spir- crbita Is  in  the  polyatomic  system 
except  these  en  the  same  atomic  center  (A).  The  effective 
one-electron  cperatcr  is  given  trore  explicitly  by 


E      (B)  E   E   -1   E 

J  =  5   q  <<»  lr   I*  > 

1   i    i  i  12  i 


(IV-2) 


where  the  q* s  are  spin-crbital  occupation  numbers  and  the 
subscript  1  is  usee!  (but  dropped  from  her<=>  en)  to  emphasize 
the  functional  dependence  of  the  operator  J  en  the 
coordinates  cf  electron  1.  When  the  eigenf unctions  are 
expanded  in  a  basis,  one  generates  the  multitude  of 
multicenter  two-electron  integrals  cf  the  type 


E   -1     E 
<  u  u  |  r   I  u  u  > . 

km   12   In 


(IV-3) 


It.  is  the  intricate  dependence  of  J  en  the  coordinates  of 
electron  1  that  does  net  allow  fcr  a  simpler  treatment  of 
this  operator,  and  so  we  seek,  an  analytic  function  with  a 
sicrpler  functional  dependence  tc  approximate  the  operator  J. 


Let  us  restrict  ourselves  tc  the  case  of  a  diatomic 
molecule  for  the  mement,  where  J  is  centered  en  a  single 
atom  (D).  Then  fcr  large  interatomic  separation  J 
approaches  the  approximate  form 


E     E  -1 
J   =  K  r 

1D 


(IV-4) 
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where  the  total  electronic  charge  is  reduced  to  a  point 
charqe  en  the  distant  atomic  center.  This  corresponds  tc 
replacinc  the  electronic  coordinates  of  electron  2  in  (IV-2) 
by  the  nuclear  coordinates  of  atom  P,  that  is, 


E      (R)  E   E   -1   E 
J   =  £    q  <<$  jr   It  >. 
i      1   1   1B   1 


(IV-5) 


At  lesser  interatomic  separation,  where  the  charge  density 
can  no  loncer  be  considered  a  point  charge,  the 
interelectronic  distance  must  undergo  an  effective  dilation 
to  account  for  the  mere  diffuse  charge.  Furthermore,  as 
electron  1  approaches  the  nucleus  of  atom  R,  the  nuclear 
attraction  should  net  witness  any  appreciable  screening,  and 
therefore  any  model  function  should  go  to  zero  as  this 
distance  goes  to  zero.  The  following  model  potential 
function,  which  we  now  choose,  has  these  desired  properties: 


_E     E     D  -»     -1 
J   =  B  T  1-J  <r   )  lr 

1B    1P 


(IV-6) 


where  X  is  a  screening  function  with  the  asymptotic  behavior 


X  (0)  -  1 


(IV-7) 


and 


X  (oo)  =  0.  (IV-8) 

Such  a  screening  function  is  inherent  to  the  Thomas-Fermi 
model  (Thomas,  1928;  Fermi,  1928)  of  the  atom.  A 
significant  feature  of  the  Thoiras-Fermi  screening  function 
is  that  it  is  universal  for  all  neutral  atoms  with  respect 
to  the  d i dens icn  less  variable 
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2  E  -1/3 
x    =  r   f  <3n/8)  Z  ] 

1E     1B 


(IV-9) 


in  terms  cf  the  nuclear  charge.  The  distance  r  is  assumed 
to  te  given  in  Eohrs,  Potentials  of  the  form  (IV-6)  are 
certainlv  nothing  new  and  most  such  potentials  (Helliran, 
1935)  aoopt  ar  exponential  tehavior  for  the  screening 
function  X.   fce  irake  the  choice 


E  -»  E  B   E  E  B  E  -» 

X  (r   )  =  ^  A  X  in    1  hi  a  ;  r   ) 

1 B     t   t  t   t  t  t  t   ID 


(IV-10) 


where 


B  — >       r.  - 1 

X  (n,l,m,a;r   )  =  r    exp  (-ax   )1   (3   ,ps   )       (IV-11) 

t  1F      1F         1E   Ira   1B   1B 

and  where  the  exponential  factors  are,  in  principle, 
universal  and  could  he  used  for  all  atoms  (with  the 
appropriate  chanqe  in  x)  .  Av.y  angular  distortions  are  taken 
into  account  with  the  real  normalized  spherical  harmonics. 
The  linear  coefficient  are  to  be  determined  by  fitting  the 
approximate  form  of  J  in  (IV-6)  to  its  true  form  in  (IV-5) 
alcno  with  the  constraint  in  (TV-7).  We  have  merely 
expanded  the  operator  J  in  a  basis  of  Slater-type  functions 
(with  a  somewhat  universal  choice  cf  exponents) . 


The  lirear  coefficients  are  computationally  quite  simple 
to  handle  since  they  can  be  taken  outside  the  integrals  in 
which  they  appear.  The  coefficients  can  te  determined  by 
substituting  (IV-6)    fcr  (IV-5)   in  the   one-center  Coulomb 
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inteqrals  cccurrinq   on  atcm   B  arid   matching  their   values, 
that  is. 


E  _E   B       B   B   P 

<d  j  j  14  >  =  <t>  i  j  |$  > 

i     i      i      i 


(IV-12) 


for  each  eiqenf u ret  ion  on  atom  B.   In  more  detail  we  have 

FE-»B  B       B       B    -»         -1       B 

N   r<*    in       I*    >    -  5    A    <<Ji    |X    (r      )r      I  <t>    >  ] 
i      1B      i  t      t      i      t      1E      1B      i 


(B)     E       B    B       -1       E    B 
=  5        q   <<D   *    |r      !0   0    > 

i         i     i   i     12     i   i 


(IV-13) 


fchere  the  Coulomb  integrals  en  the  righthand  side  cf  the 
equation  are  one-center  integrals  which  we  properly  compute 
and  retain.  fce  have  as  many  equations  like  (IV-13)  as  there 
are  eiqenf unctions  en  that  atom,  consequently  we  can  fit  up 
to  just  as  many  linear  coefficients.  The  number  of  terms 
that  are  reeded  in  (IV-1C),  though,  are  usually  much  less 
than  that  (due  to  its  relatively  smooth  behavior),  sc  seme 
type  of  (weiqhted)  least-squares  fit  procedure  should  be 
adequate.  One  ether  condition,  perhaps  even  more  important, 
that  can  be  satisfied  is  the  total  two-center  nuclear 
attraction  experienced  by  that  atom,  qive-n  by 


A  (B)  B   E   -1   B      A  B  — > 
Z  £    q  <0  I  r   |  <t>  >  =  Z  J  (r   )  . 

i    i  i  u  i        ea 


(IV-14) 


In  the  qeneral  case  of  a  polyatomic,  this  would  te  a  suit 
over  all  the  ether  atcms.  Thus,  in  addition  to  the 
equations  in  (IV-13),  we  can  have 
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B       -  1  B    B    -»  -1 

N   Tr        -  £    ft    X    (r      )  r       1 

AE  t       t    t       AE       AB 

|3)     B       E       -1       B 
=  5         q    <$     |r       |$    >. 

1  1       1        1A       B 


(IV-15) 


In   order    to   determira   the    linear   coefficients,      these    linear 
equations   can    t-<=    written    in   the    matrix    form 

T    A    =    G  (IV-16) 

where 


A      =    A    , 
t  t 


(TV-17) 


B       E       3    -»  -1       B 

T        =    l»    <♦    |  X    (r      )r       |(j)    >, 

it  i      t       1B      1E      i 


(IV-18) 


E    B    -»         -  1 
T  =    N    X     <r      )  r      , 

M+1,  t  t       AE       AB 


(IV-19) 


T  =    1    if    1      =    C, 

M+2,t  t 


=    0    if    1      >    C, 

t 


(IV-20) 


E       B       -  1       E 
G      =    N    <4>     |r       |<t>    > 
i  i      1E      i 

(B)     E       BE       -1       BE 
-5        q    <<*    $    |r      !♦*>, 

i         i     i  i     12     i  i 

B    -1  (B)     B       B       -1       B 

G  =    N    r      -  ^         q    <(J>    |r       I  4>    >, 

H*1  AB       i  1       j       1A       i 


(IV-21) 


(IV-22) 


and 


G  =    1. 

K+2 


(IV-23) 


la- 


in   qeneral    M+2    is    larger    than      the    number   of    terms    needed    in 

(17-10),    and    so    we    have    an    overdet^rmined    set.       If    W+2    and    K 

were   equal,    then   T    would    he    square    and 

-1 
A    =    T      G,  (IV-24) 

provided  that   T  is   non-sinqular  (no   linear  dependencies) . 

However,   in  the  qeneral  case,   one  can  write  the  normal  set 

of  equations  (IV-16)  a? 

t         t 
(T  wT)  A  =  IT  H)G  (IV-25) 

where  w  is   an  arbitrary  diaqonal  matrix   whose  elements  act 

as   woiqhts  tc   the  oriqinal   M+2   equations  (IV-16).     The 

choice  of  a  urit  matrix  would  correspond  to  the  usual  linear 

least-squares  fit   procedure;   nevertheless,   the   option  of 

qivinq  mere   importance  to   certain  conditions   such  as   the 

external  ruclear   attraction  (IV- 14, 15)    or  to   the  valence 

orbitals  ever  the  core  orbitals  can  be  useful.   Since  we  now 

have  square   matrices,   we  can  solve   for  A  in  terms   of  the 

qenerali2€d  inverse  of  T,  that  is, 


t    -1 
A  =  (I  wT)   (Tw)  G. 


(IV-26) 


The  exponents  in  (TV-10)  could  be  determined  once  and  for 
all  by  a  mere  fit  to  the  universal  Thomas-Feriti  screening 
function  nenticned  earlier,  but  this  would  build  in  the 
deficiencies  of  the  Thomas-Fermi  model  (particularly  for 
small  atoms).  S  better  approach,  perhaps,  was  suggested  by 
scire      work      carried      cut      by      Csavirszky       (1968,1973).  He 
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performed  a  variational  calculation  tc  determine  an 
analytical  solution  tc  the  Thomas-Fermi  equations  with  the 
following  trial  function  for  the  screening  function: 


-  3  '•  V   ? 


(IV-27) 


X(x)  =  (A'e     +  A"e     )  . 
His  calculation  yielded  the  exponents 

a1  =  C.178255C   and 

a"  =  1.759339  (IV-28) 

as  well  as   values  for  A'  and   A"  {tut  these  are   net  cf  any 

concern  tc  us).   This  corresponds  tc  the  following  choice  of 

exponents  in  cur  model  potential  expression  (TV-10): 

a   =  2a'  =  0. 35651  18, 

1 
a   =  ai  +  a"  =  1.937591,   and 


a   =  2a"  =  3. 518678. 
3 


(IV-29) 


These  valtes  »pr?  derived  for  spherically  symmetric  neutral 
atoms  in  *he  Thomas- Fermi  model  with  somewhat  modified 
boundary  conditions  iirplied  by  the  particular  choice  of 
(IV-27) .  In  order  to  allow  for  corrections  to  the  Thomas- 
Fermi  model,  radial  and  angular  distortions,  and  deviations 
frcm  neutrality  that  would  exist  for  an  "atom-in-a- 
molecule",  some  additional  flexibility  (other  than  the 
variation  of  the  linear  coefficients)  might  be  necessary. 
Augmenting  expression  (IV-10)  with  an  extra  term  or  two 
would  protably  h«  sufficient  for  most  cases,  but  angular 
distortions  would  certainly  necessitate  the  inclusion  of 
seme  polarization  +erms.     Ore  could  also  just   choose  seme 
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qeceral  veil- rounded,  flexible  basis,  such  as  ar.  even- 
tempered  (Bardo  and  ?,  uedenberg,  19*73)  set  of  exponents,  with 
the  hopes  that  the  use  of  the  diirensionless  variable  x  and 
well-determined  linear  ccefficierts  would  he  sufficient. 


CHAPTEE  V 

AN  AEFECXIMAT*  NCNLCCAL  TNTEPATGMIC 
EXCHANGE  POTENTIAL 


In  the  last  chapter,   we  substantially  reduced  the  number 

of  two-center   two-electron  Coulomb  integrals   by  exploiting 

the   localized  picture.     We  shall   new   attempt  a   similar 

reduction  in  the  analoqous   exchange  integrals.  Again,   we 
wish  to  approximate  the  effective  one-electron  operator,  hut 

this  time  the   operator  is  ncnlccal  and   requires  a  somewhat 
different  approach. 

The  twe-certer  excharge  integrals  are  of  the  form 


<4  IK  14  > 

i   1   i 


(V-1) 


where  K  is  the  effective  nonlocal  one-electrcn  operator 
(subscript  1  dropped  from  here  en)  for  the  exchange 
interaction  with  all  the  other  spin-orbi  tals  {of  saire  spin) 
in  the  pclyatcmic  system  except  these  associated  with  this 
center  (A).   The  operator  K  is  defined  such  that 


E  (B)  B   E   -1       B 

K  *  (1)  =  5    q  <♦  |r   I*  >  *  (1)  . 

i      i     i   i   12   i    -j 


(V-2) 
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Let  us  consider  the  asymptotic  behaviour  of  (V-1)  when  the 
averaqe  i rt erelectrcnic  distance  approaches  the  internuclear 
distance,    that    is,    v. hen 


E  -1        (P)     B       D  B 

<<J    |K     |0    >    =    r        ^         q    <$    |«>    ><♦     !$    > 
i  i  API  1      j      i        i      1 


(V-3) 


E     -1    (B)  B   B    B 
K   =  r    £    q  Id  ><*  I  . 

ab  i    i  i   i 


<V-4) 


Since  K  is  nonlocal,  we  can  no  lonqer  use  any  of  its  "local" 
properties,  such  as  its  expectation  values  with  respect  to 
the  eiqe Ef unctions  en  its  own  center,  to  help  determine  an 
intermediate  fortr  fcr  its  approximation.  Cur  endeavor  thus 
far  has  leen  to  retain  a  theory  which  is  valid  throuqh  at 
least  secend-crder  in  diatcmic  overlap,  and  so  it  is  felt 
that  retaininq  K  in  its  asymptotic  form  (V-U)  ,  which  is 
second-order  in  overlap,  would  net  be  inconsistent  with  any 
ether  ap r rcximat ions  made  thus  far.  In  such  a  case,  the 
interatomic  exchanqe  can  be  treated  with  the  same  ease  as 
the  pseudepoten t ia 1  in  chapter  III. 


This  approximation  may  appear  to  be  a  bit  crude  at  first, 
but  the  naqnitude  of  these  interactions  in  this  localized 
schetre  is  relatively  small  due  tc  the  nearly  electrostatic 
nature  of  the  model.  It  should  be  emphasized  that  the 
approximations  used  fcr  the  interatomic  two-electron 
interactions  could  only  be  realized  in  such  a  localized 
separable  picture. 


CHAPTER    VI 

A    KOEEL    HAKIITONIAN     AND     AN    APPROXIMATE 
ELECTRON     ERCPAGATOF 


According  to  the  Heisenberg  equation  of  motion  fcr  the 
electron  propagator  (Lirderberg  and  Chrn,  1973),  in  the 
eneray   representation,    we    have    in    matrix    form    that 


G      (E)    =   <<a    ;a    >> 

ii  i      1      E 


-1  +  ♦ 

=    E      r  <ra    ,a    ]    >    +    <<[a    ,H  ];a    >>     ] 
i      1    ♦  i  j      E 


(vi-1) 


where 


a     •=    <<J    |qi      >      and 

i  i      cp 


a      =    <l|       I  4    > 

i  op      i 


(VI-2) 


are   the    annihilation   and      creation   operators,      respectively, 
i.ith    the    anticorcmutation    relations 

+      + 

Ta,a]      -   [a    ,a    ]      =0      and 
i      i    ♦  i      1   + 


[a    ,a    1     =  [  a    ,a    1 
i      1   ♦  i      1   ■ 


s        =    <♦    !♦   >, 

ii  i      i 


(▼1-3) 
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defined  in  terms  of  field  operators  and  a  spin-crbital 
basis.  The  irany-  e  lect  ron  Ha  iri  ltcrian,  in  second  quantized 
fciir,    is 


+  +    f 

H    =  £    h      aa      *    1/2  £         (kl|  i  j)  a    a    a    a 
il    i  1    i    1  i jkl  i    k    1    -j 


where 


h         -    <«    | k     |4    > 

il  i       1      i 


an  d 


-1 
(kl|ii)    =    <<J»    4    |r       U    *   > 

k    i       12      1    i 


-1  -2 

G       (E)     =    E       S  +    F       F  +... 

i  1  i  1  i  j 


(VI-4) 


(VT-5) 


(VI-6) 


ar«    the      usual    ere-      and    two-electrcn      integrals.         Equation 
(VI-1)    can    be    iterated    to    yield 


(VI-7) 


where 


il 


<rra    ,H],a     1    > 

i  1    ♦ 


(VI-0) 


is    the    first    ncirent   and    so    en.       Substituting    *ho    Hamiltonian 
(VI-4)    into    (VI-8)  ,    the    first    mement    becomes 


F        =    h        ♦  i   r  (ij|kl)-(il|k  j)  ]<a    a    > 

ii  ij        kl  k    1 


(VI-9) 


which  has  the  same  fcrm  as  the  effective  cne-electron 
Hatriltonian  as  originally  presented  by  Fock  (1932)  and  will 
henceforth  be  called  the  Feck  matrix.  This  more  general 
approach  to  Hartree-Fcck  theory  will  be  taken  up  in  the  next 
chapter. 
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Suppose  now  that  irstead  cf  using  the  correct  many- 
electron  Fatr.i  It  cnian,  one  substitutes  into  eguation  (VI-8) 
an  approximate  one-electron  model  Hamiltonian  of  the  fern 


H  =  £  £  a  a   =5Rn 

k   k  k  k    k   k  k 


(VI-10) 


where  the  n's  are  -just  the  occupation  number  operators  and 
the  S's  are  real,  negative  energies  of  the  noninteract ing 
electrons.  Ir  cur  present  treatment,  we  would  choose  as  our 
spin-crbital  basis  the  nonort hogonal  atomic  spin-or ritals 
generated  froir  each  projected  Hairi  Itcrian  calculation,  and 
fcr  the  energies  we  would  choose  the  ecrrespend ing 
eigenvalues,  with  pcssible  modification  (to  the  positive 
eiqenvalues,  for  example).  The  Fock  matrix  then  takes  the 
approximate  fcrir 


F    =  <rr  a  ,H],a  1  > 
ii      i      i  + 


=  5  s  e  s 

k   ik  k  ki 
or  the  matrix  form 


(VI-11) 


F    =    SES 

where   6    is      a    diagonal    matrix. 

form    tecciies 

-1  -2 

G(l)     -    E       S     *    E       SES    +    .     .     . 

-1  -  1 

=     (ES       -6) 


(VI-12) 
Eguation       (VI-7)       in    iratrix 


(VI-13) 
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where  one  car  see  that  the  poles  cf  the  propagator  will 
occur  at  the  eiqenvalues  of  the  Fock  matrix,  that  is,  at  the 
zeros  of  the  secular  determinant 

iSSS-FS|.  (VI-14) 

The  ccrrespondin c  eigenvalue  prcblen 


FC  =  SCE 
can  te  reduced  tc  the  simpler  one 
t*C    =  C'E 

where 

4  + 

C»  C  =  C'C       =1 


(VT-15) 


(VT-16) 


(VI-17) 


F»  =  -XSX 

with  X  being  the  diagonal  matrix  with  elements 

1/2 
X         -     (-6    ) 
kk  k 


(VI-18) 


(VI-19) 


This  can  be  considered  as  the  di agonalization  of  a  Feck 
matrix  in  a  basis  which  is  energy- weighted  Iflwdin 
crthoqonalized,  and  this  is  the  reason  for  the  term  "Energy- 
Weiqhted  Paxiroum  Overlap"  (F.WMO)  used  ty  Linderherg  et  al. 
(1976)         to    describe      this    method.  The    molecular      orbital 

coefficierts    in    the    original    basis    are    given    by 
V2 


C         =    C»     (E    /S    ) 

ki  ki      i      k 


(VI- 20) 

The  electron  propagator  in  this  orthogonal  basis  becomes 

-1        -1 
G(E)  =  X   (El  +  XSX)   XS,  (VI-21) 
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and  the  elements  cf  the  corresponding  one-electron  reduced 
density  tratrix  are  given  by  the  contour  integral  (Linderberg 
ft  al.,  1976) 

-1, 


<a 


-v 

a  >  =  (2ni)   )  d?  <<a  ;a  >> 


k  1 


k   1   F 


1 


=  (2-ni)   )  dF  G  (?) 
C     kl 

-1(       -1         -1 
=  (2ni)   )  dF  (X   (F1*XSX)   XS) 

C  kl 


(VI-22) 


where  C  encloses  the  occupied  energies.  An  appropriate 
defiriticr  of  the  charge  and  bond-order  matrix  *culd  in  this 
case  be 


-  1 


P  =  (2 


ni)   )  dF  |EU1 


-1 


SX) 


(VI-23) 


This  choice  leads  to 


q   =  E 
k     kk 


5  <n  >|C«  | 

i    i    ki 


(VI-24) 


for  the  formal  spin-orbital  charges  in  terms  of  the 
molecular  orbital  occupations.  This  definition  is 
equivalent  to  the  orbital  populations  according  to  Kulliken 
(1955)  . 


In  the  case  cf  positive  energies  in  (VT-10),  one  would 
have  tc  resort  to  some  other  type  of  orthogonal 
transformation  in  order  to  avoid  a  complex  Fock  matrix 
(VT-18).     Lfiwdin's   (1970)     symmetric   orthogonalisation 
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procedure  would  most  closely  resenble  the  procedure  just 
described.  The  resultinq  population  analysis  is,  of  ccurse, 
dependent  on  the  type  of  crthcgcnal  transformation  used 
unless  the  analysis  is  done  with  respect  to  the  original 
nonorthoqcnal  basis.  Population  analysis  ir  a  nonorthogonal 
basis,  however,  can  lead  to  unreasonable  charges  (negative 
or  greater  than  one).  One  could,  of  course,  always  solve 
the  complex  eigenvalue  problem  and  still  use  the  atcve 
procedure,  but,  in  general,  the  irclusicn  of  positive-energy 
(virtual)  spir-orbi tals  is  not  a  straightforward  task. 


CHAPTER  VII 
THE  GBJNfl-CANCNICAl  HARTREE-FOCK  SCHEME 


The  notion  of  an  atom  in  a  molecule  led  Slater  (1970)  tc 
consider  an  en<?rqy  functional  obtained  as  an  average  over 
multiples  arisinq  from  a  given  configuration.  His 
extension  of  this  idea  tc  configurations  with  fractional 
occupations,  known  as  the  Hy per-Hartree-Fock  method, 
however,  has  been  subject  to  some  criticism  because  cf  the 
appearance  of  cff-diaqcnal  Lagrangian  multipliers  and  lack 
of  certain  conceptual  grounds  (Jjargensen  and  flhrn,  1973)  . 
In  this  same  spirit  we  wish  to  introduce  an  average  based  on 
the  statistical  mechanical  concept  of  an  ensemble  (Abdulnur 
et  al.,  1972).  The  particular  choice  of  a  Grand  Canonical 
ensemble,  as  defined  by  its  density  operator,  leads  to  the 
elimination  of  off-diaqcral  Lagrangiar  multipliers. 


A   system  cf   nonin teract ing  electrons   described  by   the 
Hairiltcnian 


H  =  $  E  n 
i   i  i 


(vn-1) 
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can  te  described   by  the  Grand  Canonical   partition  function 
related  tc  the  density  operator 


p  =   tt  (1-n  +2  n  )/  (Hz  ) 
i      i   i  i       i 


where 


z   =  expf-  <&  -u)  /T  ]■ 
i  i 


(VTI-2) 


(VII-3) 


The  parameters  /i  and  T  are  the  thermodynamic  chemical 
potential  and  absolute  temperature.  The  average  value  of  an 
operator  A  fou  such  a  system  is  then  given  in  terms  cf  its 
trace  with  respect  to  the  density  operator,  that  is, 

<A>  =  IrfAp}.  (vII-4) 

In  particular,  the  average  value  cf  the  number  operator  is 


<n  >  =  z  /  |1  +  z  )  . 

i     i      i 


(VII-5) 


Usinq  this  ir  expression  (VTI-2)  gives  a  unique  way  of 
defining  the  density  operator  from  a  given  set  of  occupation 
nuirters, 


p   =   n    r  1-<n    >+(2<n    >-1)n     ], 
i  i  i  i 


(VII-6) 


without  any  reference  to  any  thermodynamic  parameters. 
Althouqh  there  is  nc  connection  with  eguilibrium  situations 
in  statistical  mechanics,  we  shall  refer  to  this  as  a  Grand- 
Canonical  density  operator,  which  acts  as  a  formal  device 
for  the  formation  of  averages  only.  Except  for  the  case  of 
all  inteqral  occupation  numbers,  this  form  of  the  density 
operator  qives  a  nonzero  width  (uncertainty)  for  the  total 
number  operator. 
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If  cne  considers  an  electron-conserving  (Cancr.ical) 
averaqe  value  cf  the  na ny-e lee tron  Hamiltonian  in  the 
crthcnormal  basis  which  diagonalizes  the  density  matrix,  cne 
has 


<H>  =  ^  h  <n  >  +  1/2  5  (J   *K   ) <n  n  > 

iii        ijijij    i  i 


(VIT-7) 


where  the  only  thirg  which  is  not  explicitly  determined  is 
the  averaqe  cf  the  number  operator  product  in  the  secend 
term.  If,  however,  the  averaqe  is  performed  with  respect  to 
the  Grand-Cancnical  density  operator  then 

<n  n  >  =  <r  ><n  >,  (VTI-8) 

i  i      i   i 

and   equation    (VII-3)    can    be    written 

<H>    =  ^  6    <n    >  -    1/2  ^    V       <n    ><n    >  (VII-9) 

iii  iiiji         j 

where   the   effective   one-electron   energy   is   given   by 


5      =    h      +   5    y      <n> 
i  i         1      ij       i 

and  the  effective  interaction  enerqy  by 

V    =  J   -K   . 
ii    ii  ii 


(VII-10) 


(VII-11) 


Cne  reccqeizes  that  these  energies  can  be  interpreted  as  the 
first  anc  secend  partial  derivatives  cf  the  total  energy 
functional  <H>  with  respect  to  the  occupation  numbers  as 
follows: 


a<H>/d<n  >  =  E 

i      i 


(VII-12) 


and 
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3  <H>/d<n  >d<r  >  =  3e  /a<n  >  =  V   . 
i     i      i     i     ij 


(VTI-13) 


Such    interpretations    lead    to    the    finite-difference 
approximations  to  ionization  energies 


(i) 

A  E  =  S  A   <n  > 

i     i 

and  excitation  energies 

(i->i) 


(VII-14) 


A  E 


5^<n>    +    5^J<n>    +    V      Z\  <n    M  <T1    >(VII-15) 

i     i      1     J      ij     i      1 


for  which  Kocpmans'  (1935)  theorem  is  a  special  case. 
Application  of  the  mean  value  theorem  suggests  that  the  test 
approximations  to  these  quantities  would  be  obtained  by 
determining  the  S*s  and  V's  at  some  intermediate  point,  such 
as  <n>=1/2.  Such  a  scheme  has  been  termed  the  transition 
state  method  (Slater  and  Wood,  1971)  when  applied  to  similar 
enerqy  furcticnals. 


The  effective  one-electron  (Grand-Canonical)  operator 
havinq  ttese  cne-electrcn  energies  as  its  eigenvalues,  has 
matrix  elements,  in  an  arbitrary  nonorthoqona 1  basis,  cf  the 

form 

f 
I    =  h    +  5    r  (  k  1 1  k'l')  -(kl»  |k'l)  ]<a   a   >     (VTI-16) 
kl    kl    k'l1  kf  1* 

where  the  density  matrix  elements  are  given  by 


<aa>=^c   <n>C 
k  1     i   li   i   ik 


(VII-17) 
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and  C  is  the  unitary  transfer iraticn  matrix  in  this  basis 
which  diaqonalizes  F,  ir.  ether  words,  its  columns  are  the 
associat.se  eigenvectors  cf  P.  The  iterative  self -consistent 
construction  and  solution  of  equation  (VII-16)  is  what  has 
come  to  be  called  the  Grand-Canonical  Hartree-Fock  (GCHF) 
method,  viith  the  total  energy  given  by  equations  (VII-9), 
(VII-10)  ,  and  (VII-  11)  . 


CHAPTER  VIII 

B  SELF-CONSISTENT  CHARGE  ANC 
CCNFIGURATTCN  PRCCEDORF 


Up  tc  now  we  have  said  very  little  about  the  manner  in 
which  the  spin-orbital  occupations  are  assigned.  This  is 
scirewhat  cf  a  sensitive  subject,  here  arc!  in  many  other 
"bcildinq-blcck"  approaches.  This  question  has  its  roots  in 
the  early  days  of  valence-bond  methods.1  We  would  like  tc 
make  an  analcqy  tetween  our  trethcd  and  the  single- 
confiquration  cf  nonorthoqona 1  crbitals  method.1  This 
method  amounts  tc  a  condensation  cf  many  configurations  to 
one,  built  up  from  distorted  (hybridized)  atomic  orbitals, 
which,  in  essence,  is  wha«-  Coulscn  and  Fischer  (1949)  did 
for  the  rydroqen  irclecule.  The  application  of  this  method 
to  larqer  systems,  however,  is  quite  limited  because  cf  the 
difficulty  in  calculatinq  the  iratrix  elements  of  the 
Hairiltonian,   sir.ce   there  is   no  orthogonality   between  the 


*The  literature  on  these  methods  is  -just  too  extensive  tc 
review  here,  tut  we  sucqes*  Gerratt  (1974)  for  an  excellent 
account  cf  the  subject  as  it  applies  here. 
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orbitals.  The  model  itself  provides  a  sensible 
interpretation  of  chemical  bondinq  in  terms  of  a  distortion 
cf  the  participat inq  atomic  orti+als  combined  with  a 
recouplirc  of  the  spins.  The  method  has  been  shown  to  be 
capable  cf  yielding  useful  estimates  of  molecular  binding 
enerqies  and  cthpr  properties.  foe  view  our  method  as  a 
direct  ireans  cf  cbtaininq  such  distorted  ncnort hcgcra 1 
orbitals,  wtere  the  orthcgcna  lit  y  problem  has  been 
incorporated  into  the  pseudopotentials.  The  renor realized 
pseudoor t itals  (nncenst  rained  furc^ons  cf  chapter  IT)  then 
act  as  these  hybrid  orbitals.  The  tctal  enerqy  we  cempute 
frcm  these  orlitals  is  valid  thrcuqh  second-order  in  overlap 
(from  the  flulliken  approximation)  and  so  we  should,  within 
the  int^qral  approximations  made,  be  able  to  irimic  such  a 
scheme,  es  far  as  we  can  tell.  There  at?  some  other  basic 
differences,  however.  These  methods  generally  separate  the 
space  and  spir  parts  of  their  wa vef unctions ,  and 
consequently  couple  their  spins  with  proper  generalized  spin 
functions. 


The  cencept  of  "atciric  valency"  arises  when  molecules  are 
allowed  tc  separate  into  their  constituent  atoms.  When  the 
most  qeneral  linear  cemtiratien  of  spin  couplings  is  formed, 
the  resultinq  separated  atomic  states  can  be  regarded  as 
well-defined  "valence  states".  Any  less  ccuplirg,  however, 
can  lead  tc  ill-defined  r.cnccupled  states,  particularly  when 


-52- 


degenerate  state?  ate  involved.  This  we  feel  is  mostly  due 
to  the  separation  of  space  and  spin,  connected  with  an 
imbalance  of  synnetry  constraints  cr  the  two  parts.  In  our 
treatment  we  use  spin-orbitals ,  which  frcm  *he  beginning 
puts  space  and  spin  en  similar  footing.  A  single 
configuration  of  spin-orbitals  will,  in  general,  not 
correspond  to  a  proper  spin-coupled  state,  but  the  atcnic 
"valence  states"  remain  rather  well-defined,  since  each 
valence  electron  is  accomodated  in  a  distinct  (apart  from 
exchange)  atonic  orbital  with  a  distinct  spin. 


One  desires  a  configuration  which  allows  one  to  describe 
the  disscciation  process  as  a  smooth  recoupling  of  the 
orbitals  from  a  per  feet ly- paired  state  to  the  proper  atomic 
states.  Particularly  attractive,  in  this  regard,  is  the 
spin-valence  theory  employed  by  Heitler  (1934),  where 
eiqenf unctions  of  the  atomic  Hamiltonians  are  ccupled 
toaether  to  fcrm  the  molecular  state.  Tn  fact,  this  was  the 
motivation  behind  Mcffitt's  (1951)  atom-in-molecule 
approach.  The  basic  question  is  "hew  do  we  choose  cur  spin- 
orbital  cccupaticns  in  our  single-configuration,  which  has 
these  desired  properties,  without  falling  into  the  traps 
that  ether  methods  do  when  all  the  proper  couplings  ate  not 
taken  intc  account?"  Before  we  answer  this,  let  us  mention 
two  methods  which  come  close  in  spirit  to  what  we  are  abcut 
to  describe.    The   first  is  Hund's  rule   coupling  (Getratt, 
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1971),  and  the  second  is  alternate  molecular  orbital  (AFO) 
theory  (Faunc2,  1967),  We  mention  them  because  of  the 
physical  picture  they  present;  however,  in  th^se  two  related 
methods,  space  and  spin  are  separated,  and  the  resultant 
spin  functions  are  complicated.  We  feel  that  if  one  instead 
uses  spir-or ri tals  and  relaxes  seme  of  the  constraints1 
normally  imposed  on  the  spin  (Gunnarsson  et  al. .  1977),  then 
the  same  physical  picture  suggests  a  very  simple  scheire  for 
makinq  our  assiqrmerts. 

The  procedure  is  to  start  out  with  a  configuration  which 
corresponds  tc  a  covalent  structure  that  one  generates  from 
valence  shell  e  lectrcn-pair  repulsion  theory  (VSFIPRT)  as 
described  in  any  general  chemistry  text.  The  spins  are 
assiqned  such  that  the  electrons  on  each  atom  are  ccucled  tc 
qive  maxinuni  resultant  and  such  that  the  resultant  spin  on 
each  atoo  alternates  in  sign  with  respect  to  each  cf  its 
bonded  neiqhbcrs.  This  should  be  done  such  that  the  total 
resultant  spin  is  the  desired  cne  (if  possible)  .  Consider 
the  example  of  a  carben  menoxide  irclecule.  It  has  the  valid 
structure 

i  c  si  c  (vii i-1) 

3D. 

correspendinq  tc  the  co n f iqurat icn 


'These  constraints  were  necessary  in  the  other  treatments 
to  qet  prcper  separation  of  the  atomic  states.  This  is  net 
necessary  in  cur  sitrple  treatment. 
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(C1s)  (C1s«)  (C1s)  (els')  (02s)  (C2s')  (C2s)  (C2s') 

(C2px)  <02px')  |02pz)  (C2pz«)  (02py)  (C2cy')  (VIII-2) 

where  the  priires  denote  opposite  spin.    There  are  cf  course 

ionic  structures   as  well,    but  we   shall  incorporate   them 

later.   The  main  point  here  is  that  this  structure  separates 

into  the  atomic  states  represented  by  the  configurations 

(C1s)  (C1s»)  (C2s)  (02s')  (0?px)  (02px')  (02py)  (02pz)  (VITI-3) 
and 

(C1s)  (C1s'}  (C2s)  (C2s»)  (C2pz')  (C2py*)  (VIII-4) 

each  satisfyirq  Hund's  rule.  This  type  of  structure  can  be 
qeneratec  for  any  systeir  where  the  atcms  can  be  divided  into 
the  subsets  in  such  a  way  that  no  two  atoms  which  belong  tc 
the  same  subset  are  neighbors  tc  each  other.  Such  a  system 
is  called  an  "alternate"  system.  For  nonaltemate  systems 
the  situation  can  be  quite  different.  In  fact,  in 
nonaltemate  solids  this  can  give  rise  tc  anti ferromagnetic 
conduction  sheets.1  Tn  cyclic  molecules  (the  only  molecules 
which  car  be  nonal ternate)  ncnalternacy  is  often  accompanied 
by  unusually  stable  ions,  radicals,  "sandwich"  complexes,  or 
some  higher  level  cf  alternacy.2  Even  though  this  would  be 
an  interesting  topic  on  its  own,  we  shall  not  dwell  on  it 
here.  The  main  point  we  wished  to  make  here  is  that  withir 
each  atoii  we  have  parallel   alignirert  cf  spins  (Hund's  rule, 


*Fcr  example,  an  FCC  lattice  might  have  alternating  planes 
cf  parallel  spin,  thus  a  higher  order  of  alternacy. 
2Such  as  sicma  and  pi  planes. 
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f e rrcmaq retic  couplinq)  and  between  atoms  we  have  (at  least 
for  alternate  systems)  antiparallel  alignment  of  spins 
(antif errcmagne t ic l  couplinq).  This  is,  of  course,  just  a 
qeneral  cuidelire,  art!  one  can  make  any  occupation 
assignments  ore  wants,  provided  that  care  is  taken  not  to 
break  symt^trv  with  exchanqe  polarization  (the  scheme  just 
describee!  will  net  break  symmetry  vhen  applied  to  alternate 
systems)  . 

Sc  far  we  have  just  devised  a  way  to  assign  spins  such 
that  we  cet  strccth  uncoupling  upcr  separation.  Suppose, 
thouqh,  that  this  sinq]e  neutral  structure  does  not  prcperly 
describe  the  molecular  state  or  the  separated- atom  states. 
What  if  ve  krew  that  the  irclecule  had  a  triplet2  qrcund 
state  which  can  not  be  formed  by  a  smooth3  coupling  cf  the 
separated-atorr  qrourd  triplet2  states?  what  if  we  suspect 
the  qrourd  state  to  be  a  qreat  deal  ionic  in  character? 
These  types  cf  questions  involve  what  is  perhaps  best 
described  as  the  crcssinq  of  different  single-configuration 
diatatic  states  (O'lalley,  1971).  The  simplest  prodedure 
would  be  to  compute  both  of  the  diabatic  states  involved  and 
"uncross"  them.   In  fact,  for  the  first  case  where  the  curve 


1Nct  tc  exclude  the  possibility  cf  complete  ferromagnetic 
ccuplinq. 

2Ey  "triplet"  we  mean  two  upaired  parallel  spins;  total 
spin  is  never  properly  taken  into  account  in  our  treatment. 

3"smootb"  meaninq  nc  sudden  charge  transfer  or  spin  flips. 
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crossing  corresponds  to  a  spin  flip,  this  is  atcut  the  crly 
alternative  ore  has  in  the  present  treatment;  however,  in 
the  second  case  we  only  hav?  charge  transfer  (no  spin  flips) 
takinq  place.  This  brings  us  to  the  main  theme  of  this 
chapter  -  a  self-consistent  charge  and  configuration 
procedure . 


As  mentioned  earlier  there  are  other  valid  (ionic) 
structures  one  can  write  down  for  carton  monoxide  besides 
the  covalent  cne  in  (VITI-1).  These  ether  structures  are 
nothing  irore  than  charge  transfers  frcm  the  covalent 
structure  (dismissitq  spin  flips).  By  employing  fractional 
occupations  ir  cur  configuration  we  can  go  smoothly  frcm  one 
structure  to  another  (one  diabatic  state  to  another)  .  In 
aeneral,  the  molecular  state  itself  right  best  be  described 
by  such  a  configuration  with  fractional  occupations 
(configuration  mixinq)  which  goes  "adiabaticall yM  into  the 
separated  atonic  spates  with  integral  occupations.  Hew  dees 
one  determine  these  fractional  occupations?  Sc  far  we  have 
not  even  said  what  determines  the  irclecular  state  in  our 
treatment.  Cne  could  just  define  it  in  +  erms  cf  a  given 
(valence- bond)  structure  constructed  from  the  atomic 
crbitals,  but  most  molecules  probably  can  not  bp  well- 
represented  by  a  sinqle  structure.  In  chapter  VI  we 
developed  a  means  cf  generating  mclecular  (delocali zed) 
orbitals  from  a   model  Harail^onian  built  up   froir  the  atomic 
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ortitals  and  energies.  Perhaps  the  best  way  to  describe  the 
molecular  state  then  is  in  terms  of  these  molecular 
orfcitals.  We  could  cheese  integral  occupations  for  these 
molecular  orbitals  corresponding  tc  seme  desired  state  and 
then  let  the  population  analysis  determine  the  fractional 
cccupaticns  fcr  the  atomic  orfcitals  which  gc  intc  our  single 
valence- tend  structure.  Such  a  procedure  would  te  very  much 
like  the  hype r-Hartree-rock  method  proposed  by  Slater  et  al. 
(1S6?)  fcr  crystals.  This  type  of  procedure  would  allcw  one 
to  compute  a  potential  energy  surface  which  by  censtruction 
corresponds  tc  a  molecular  state  with  a  specific  spin1  and 
angular  ircmentum1,  as  well  as  to  a  specified  set  of 
separ ated-atom  states  with  the  sane  net  spin  and  angular 
momentum.  Perhaps  the  best,  way  tc  illustrate  this  is  by  an 
example,  and  the  next  chapter  should  serve  this  purpose. 


lCnly  the   compenent  alcng  scire   common  preferred   axis  is 
explicitly  and  uniquely  specified. 


CHAPIFR  IX 
APPLICATIONS 


A.  Nitric  Cxide  Molecule 


Despite  the  important  role  of  nitric  oxide  and  its 
positive  ion  in  the  upper  atmosphere  and  its  ecologically 
undesirable  presence  in  the  exhaust  emissions  of  cur  ever- 
sc-popular  automobiles  down  here  or  earth,  there  has  teen 
relatively  little  theoretical  electronic  structure  work 
reported  for  this  first-row  diatomic  molecule.  Restricted 
Hartree-Fcck  and  configuration  interaction  calculations  have 
been  carried  cut  for  the  ground  <X2n)  and  first  excited 
(A22*)  states  of  nitric  oxide  at  their  equilibrium 
geometries  by  nreen  (1972,  197?)  to  yield  cne-electron 
properties  with  rather  litrited  srccess.1  Put  tc  our 
knowledge  ther°  has  been  no  unrestricted  Hartree-Fcck 
calculation  (spin- prelected  cr  otherwise)    reported  on  this 


1See  also   the  natural   orbital  calculation   of  Kouba   and 
Ohrn  (197  1)  . 
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odd-electron  (paramagnetic),  open-shell  system  tc  assess,  on 
a  one-electron  basis,  its  complex  photoe lectrcn  spectrum 
(Turner  et  aj.,  1S7C)  .  Configuration  interaction 
calculations  for  the  positive  icn  have  been  done,  however, 
in  an  attempt  to  assiqn  the  various  observed  static  (11 
between  9  and  2H  eV)  by  Lef  ehvre-B  ticn  and  floser  (1966)  and 
later  by  Thulstrup  and  Chrn  (1972)  .  The  assignments  appear 
to  be  resolved,  althouqh  alternate  assignments  have  been 
proposed  ty  Collin  and  Natalis  (1968). 


fie  make  no  attempt  here  to  confirm  or  challenge  these 
assiqnmerts,  hut  rather,  tc  test  cur  theory  for  its 
strenqths  and  weaknesses  on  a  reasonably  small  but  complex 
system  before  applyinq  it  to  a  much  more  ambitious  problem. 
Since  our  aitr  is  tc  look  into  the  chemisor ption  of  nitric 
oxide  en  a  nickel  surface,  perhaps  we  should  know  hew  the 
theory  works  on  nitric  cxide  itself.  From  the  outset  we 
expect  (cr  at  least  hope)  that  the  dissociation  process,  in 
terms  of  the  interaction  of  localized  atomic  crbitals  and 
the  associated  total  energy,  should  be  reasonably  well 
described,  since  the  total  energy  is  not  dependent  upon  how 
the  total  density  is  broken  up  into  various  contributions 
but  rather  on  the  to+al  density  alone.  On  the  ether  hand, 
properties  associated  with  the  delocalized  molecular 
crbitals  derived  frcir  our  model  Hami  Itcnian  ,  such  as 
KccplIans,  energies,   mus*  be  taken   with  less  reliance.    Tn 
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connection  with  this,  charqe  transfer  (between  valence-bend 
structures,  if  you  wish)  is  expected  to  be  a  sensitive 
issue.  This  toe,  in  our  procedure,  depends  en  the  model 
Hamiltonian  usee  tc  generate  the  irclecular  ortitals  from 
which  the  poplulation  analysis  proceeds.  For  the  most  part, 
we  will  avoid  such  coirpl  ica  ticn  s  by  restricting  ourselves  to 
sinqle  valence-bcnd-type  structures,  that  is,  no  charge 
consistency  will  in  qeneral  be  attempted. 

Fcr  nitric  oxide  in  its  grcund  (X*n)  state,  we  will  use 
the  conf icuraticn 

(C1s)  (C1s«)  (Nls)  (N1s«)  (C2s)  (C2s»)  (N2s)  (N2s«) 

(C2p7)  (N2pz«)  <02px)  (N2px«)  (02py)  (02py»)  (N2py»)  (IX-1) 
in  terms  cf  "perfect-paired"  atomic  spin-orbita Is  where  the 
primes  denote  the  uaicrity  spin.  This  is  unitarily 
equivalent  to  a  molecular  configuration 

(16")  i-\P)  [It)   (2^«)  (3f)  (36')  (46)  (H6«) 

(5€)  (5£«)  (1nx)  (1ttx')  (1ny)  {Irry')  (2ttv«)  (IX-2) 

where  the  26",  <*6  and  2ny  orbitals  would  be  of  antibonding 
character  and  the  16  and  1-ny  ortitals  cf  rather  nonbending 
character.  In  our  procedure,  the  nitrogen  and  cxygen 
orbitals  in  (IX-1)  are  determined  separately  in  an 
alternating  fashion  ccupled  tc  each  ether  in  a  self- 
censistent  manner  as  described  in  earlier  chapters.  From 
these  atomic  ortitals,  a  model  Hamiltonian  is  constructed 
whese  natural  crbitals  (canonical  Hartree-Fock  solutions) 
correspond  roughly  to  the  molecular  orbitals  in  (IX-2)  . 
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The  basis  sets  used  for  nitrcger.  and  oxygen  are  STC's  of 
double-zeta  quality  (Roetti  and  dementi,  1<574)  plus  d-type 
polarization  functions,  for  a  total  of  15  functions  cr  each 
atom.  The  model  potentials  were  expanded  in  a  set  of  4  s- 
type  functions  and  one  p-type  function  polarized  along  the 
bond  axis  (z-direction  here).  tie  chose  an  even-tempered  set 
cf  exponents  (0. 5,  1  .C,2 .0,4. C)  scaled  by  the  cube  rcot  of 
the  nuclear  charge,  as  suggested  in  chapter  IV.  The 
exponents  for  the  tasis  sets  and  the  model  potentials  are 
listed  in  Tables  1  and  2.  No  attempt  was  made  to  optimize 
these  exponerts,  ard  the  polarization  functions  were 
selected  rather  arbitrarily. 


The  first  calculation  was  carried  cut  at  the  equilibrium 
geometry  cf  2.17  a.u.  with  the  procedure  just  described 
above  arc  in  previous  chapters,  but  the  spin-orbital 
cccupaticn  nunrbers  were  kept  fixed.  In  other  words,  any 
charqe  transfer  predicted  from  the  population  analysis 
between  each  cycle  was  ignored.  As  stated  atove,  this 
corresponds  tc  a  single  valence-bend- type  structure.  This 
is  dene  tc  eliminate  any  effects  dup  directly  or  indirectly 
*o  the  model  Harriltcnian.  The  result  was  a  total  energy  cf 
-129.44  a.u.,  which  is  about  0.10  a.u.  below  the  estimated 
Hartree-Fcck  limit  (Hcllister  and  Sinanoglu,  1966).  The 
subseguent  model  (FWMC)  Hamiltonian  calculation,  however, 
gave  molecular   crtital  energies   which  were   essentially  in 
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Table  1.   Easis  sets  for  nitrogen  and  oxygen. 
a  a 


tvp< 


nitrogen 


oxygen 


Is 

8.49597 

9.46635 

Is 

5.9eeu« 

6. £3768 

2s 

2.26  086 

2.68801 

2  s 

1 . 42457 

1.67543 

2p 

3.2  4  933 

3.69  445 

2p 

1 .49924 

b 

1.65664 

3d 

2.00000 

2.00000 

Frcn;  the  tables  of  Foetti  and  Cleicenti  (1974). 


Unoptimized  polarization  function. 


Table  2.   Model  potential  exponents  for 
nitrogen  and  oxygen. 


Type 


nitrcgen 


oxygen 


1s 

0.95647 

1.C0CCC 

1s 

1.91293 

2.00000 

1s 

3.82586 

4.COC0O 

U 

7.65173 
b 

8.COC0O 

2pz 

1.91293 

2.00000 

Even-tenpere<3  set  cf  functions  (C  .  5,  1 .  0  ,2  .  0,4  .  0) 
scaled  t v  Z**  (1/3)  . 

Polarization  function  along  bonding  axis  (1.0). 
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reverse  cf  the  ckserved  crdei  for  the  56"  and  2ny  orbitals, 
and  qave  an  ancnalcus  3^  orbital  5C  eV  below  the  observed 
value. 

In  ar  effort  to  see  what  the  source  of  these 
discrepancies  was  (assuming  that  Kcopmans*  approximation  is 
not  all  tc  blame),  a  second  calculation  was  performed  where 
charqe  transfer  was  allowed  tc  take  place,  that  is,  charge- 
consistercy  was  attained.  The  result  was  a  somewhat  higher 
total  energy,  with  some  improvement  in  the  molecular  orbital 
enerqies,  but  the  lasic  discrepancies  persisted.  The  net 
charqe  transfer  was  only  C.C15  (from  oxygen  to  nitrogen!), 
which  sugqests  that  the  single  (neutral)  valence-tend 
structure  is  perhaps  net  so  inadeguate.  Could  it  be  that 
the  model  Hacriltonian  constructed  from  just  the  occupied1 
atomic  crtitals  is  toe  peer? 

The  next  loqical  step  appeared  to  be  an  improvement  of 
the  model  Hamiltcnian.  As  mentioned  in  chapter  VTT,  the 
inclusion  of  virtual1  orbitals  in  *he  EWMO  Hamiltonian 
reguired  either  complex  arithmetic  (a  drastic  and 
undesiratle  step)  or  an  alternate  method  of  treating  the 
ncnor thoqenalit y   which  subsequently   alters  the   population 


1With  respect  to  the  model  Haniltcnian,  "occupied"  merely 
implies  a  neqative  orbital  energy  and  "virtual"  a  positive 
energy. 
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analysis.  The  second  alternative  requires  an  extra 
diagcnali zaticn  of  the  metric  matrix  (if  LOwdin's  (1970) 
symmetric  orthogcnalization  is  used)  but  is  much  easier  to 
implement  coaputat iona lly.  This  leads  to  a  population 
analysis  which  is,  in  part,  arbitrary  if  done  it  the 
orthogonal  basis  {which  is  necessary  to  quarantee  reasonable 
occupation  numbers  in  a  charqe-consistent  calculation) . 
Such  a  calculation  was  performed.  There  was  a  substantial 
improvement  tut  the  relative  crderinq  of  the  orbital 
enerqies  has  still  in  disagreement  with  the  observed  values 
as  they  beve  beer  assigned.  Charge  consistency  gave  little 
mere  improvement  t .  The  main  discrepancies  are  ir.  the 
relative  order  of  the  varicus  singlets1  and  triplets*  and 
the  apparent  reversal  cf  the  56*  and  2ny  levels. 

Cne  final  improvement,  still  involving  the  virtual 
ortitals,  was  attempted.  In  the  model  Hamiltonian,  as 
constructed  above,  the  virtual  ortitals  play  a  role 
equivalent  tc  that  of  the  occupied,  and  the  energy 
associated  with  each  virtual  orbital  should  then  be  the 
eneray  it  would  have  if  it  were  indeed  occupied.  This  is 
related  to  the  problem  encountered  in  cenf iquration 
interaction  and  many-body-perturbation   calculations,   where 


lThe  terms  "sinqlet"  and  "triplet"  refer  to  +he  final 
ionic  states  resultinq  in  the  deletion  of  a  nsajority-  or 
minority-spin  orbital,  respectively  (the  "sinqlet"  actually 
beinq  an  averaq^  of  sinqlet  and  triplet  states). 


the  Hartr€€-Fcck  virtual  crbitals  used  in  *-heir  expansions 
arise  from  a  potential  generated  by  +oo  many  electrons. 
This  erroneous  potential  comes  about  from  a  miscancella tion 
of  Cculopt  and  exchanqe  interactions  an  electron  has  with 
itself.  This  cancellation  is  an  important  artifact, 
however,  of  the  Hartree-Fock  method  which  allows  the  use  of 
a  sinqle  operator  for  all  the  occupied  {and  virtual) 
orbitals  (Huzinaqa,  et  al,  1973) .  An  approximate  first- 
crder  correction  to  this  can  be  achieved  by  just  rescaling 
the  Coulomb  and  exchanqe  contributions  to  the  orbital 
enerqies  to  account  for  the  proper  number  of  particles 
involved  in  the  potential.  For  a  charge-consistent 
calculation  this  involves  all  partially  occupied  orbitals, 
each  with  its  own  scalinq  factor.  Such  a  correction  was 
made  and  the  results  are  qiven  in  Table  3,  along  with  the 
experimentally  observed  values  and  ether  calculated  results. 
In  Figure  2,  a  simulated  photoelectrcn  spectrum  is  shewn 
where  the  profiles  are  generated  by  qaussians  of  constant 
half-width  (  FHHP!  =  0  .  ce  V)  weiqhted  by  the  spin-degeneracy  of 
the  resultant  ionic  states.  The  inset  is  a  tracing  of  the 
FSCA  result  (Sieqbarn  et  al.,  1969).  The  prevailing  feature 
is  a  shift  to  higher  energies,  particularly  for  seme  of  the 
triplet  states  relative  to  the  singlets.  In  order  to  retain 
a  total  energy  which  is  independent  of  the  model 
Hamiltonian,  the  occupations  again  were  kept  fixed  in  this 
calculation.    The  atomic  orbitals,   also  independent  of  the 
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Table    3.      Calculated    and   experiirenta  1   ionization    potentials 
for    nitric   ojcide    (X2n)  . 


Ion  a 

state        MO 


t  c 

expt        expt         e  x  p  + 


this 
calc 


calc        calc 


in 

C1s  • 

5  4  4  .  C 

564.2 

3n 

01s 

5U3.3 

565.4 

in 

N1s  ' 

tin.  e 

42c.  1 

3r 

N1s 

410.3 

422.1 

>r 

36*' 

43.8 

51.2 

'n 

3  6* 

40.6 

47.6 

in 

46*' 

23.3 

28.0 

21.0 

*v. 

46* 

21.7 

29.9 

20.0 

A*n 

56*' 

18.5 

18.3 

18.3 

18.6 

19.4 

17.8 

*r\ 

56- 

16.7 

16.5 

17.2 

16.1 

17.0 

16.9 

«s* 

1-ny' 

17.8 

18.3 

lA 

1-nx' 

19.3 

18.9 

15.8 

17.7 

17.8 

■5- 

18.4 

17.3 

17.0 

17.2 

■S- 

17.  6 

16.6 

17.0 

17. C 

M 

1nx 

16.8 

15.7 

21.9 

16.0 

16.  1 

a*S* 

1n  v 

14 

15.7 

14.0 

19.4 

15.0 

14.3 

Xl£  + 

2ny 

10 

9.  3 

9.3 

15.8 

9.3 

12. C 

Via    Koo^:Irans,    theorem. 

ESCA,    Siegbahr.    et    al.     (1969).     (Vertical) 


A 


Turner    et   al    (1970).     (Adiabatic) 
Collins    end    Natalis    (1968).     (Adiabatic) 


Lefebvre-Bricn    and    Mcser    (1966).     (Adiabatic) 


Thulstrup    anc   Obrn    ( 1 972)  .     (Adiatatic) 


Figure  2.  Calculated  and  experimental  phctoelectton 
spectra  for  ritcic  oxide  (X2-n)  .  Convolution  done  with 
qaussians  with  FfcHM =1 .5e V.  Inset  from  ESCA  (Siegfcahn,  et 
al.,  196S). 
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mcdel  Ha ailtonian ,  are  illustrated  in  Figures  3  and  4  by 
"hidden- line"  plcts  (Mcintosh,  1975)  in  the  yz-plane  for 
nitrogen  and  cxyqen  respectively.  The  inclusion  of  virtual 
orfcitals  in  th»  model  Familtonian  tends  tc  be  a  sensitive 
issue  as  can  te  the  charqe  transfer,  tut  unlike  the  charge 
transfer  problem,  it  has  no  direct  effect  en  the  total 
enerqy  or  the  atomic  orbitals,  provided  the  occupations  are 
restricted  to  zero  and  unity. 

At  this  point,  cne  might  begin  to  have  some  doubts  abcut 
the  various  apprcxi  ira  tions  used  ir  cur  treatment,  or  the  use 
of  a  single  nonsymmetric  spin-polarized  valence-bend 
configuration.  One  must  reserve  judgement,  however,  until 
each  of  these  approximations  has  teen  carefully  examined  on 
a  variety  of  systeirs,  for  one  can  not  help  but  single  out 
the  model  Hamiltcnian  as  still  a  prime  target  of  scrutiny. 
One  miqht  alsc  te  suspect  of  Kccpraans'  approximation  for 
this  system.  A  spin-projected  unrestricted  Hartree-Fock* 
calculate  en  wculd  help  in  this  regard. 

Let  us  now  leck  at  what  should  te  a  more  premising  side 
of  this  story.  The  total  energy  that  we  obtained  at  the 
equilibrium  geometry  (vith  fixed  occupations)  was  about  6  eV 


2A  spin-polarized   X-alpha  calculation  of  Connelly   et  al . 
(1973)  is  in  reasonable  aqreement  kith  experiment. 

lUsing  the  same  basis  and  Grand-Canonical  Hart ree-Fock . 


Fiqure    3.         "Hidden-line"    flots   of   the    atomic    crtitals 
on    the    nitrcqen    atcir.       yz-plane     (x=C.O). 
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Fiqure  U.         "Hidden-line"  plcts  of  the  atomic  ortitals 
en  the  os^qen  atom.   yz-plane  (x=C?). 
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lover  than  that  of  the  separated  Hartree-Fock  atoms.1  The 
experimental  dissociation  energy  is  6.5  eV.  This  is  either 
fortuitous  or  very  enccuraqinq.  A  full  potential  energy 
curve  is  certainly  in  order  at  this  point.  Fcr  the  most 
rapidly  ccnvercent  procedure,  ere  starts  the  potential 
enerqv  curve  at  the  separated-atcm  limit  where  one  is 
quaranteec  the  proper  limit  (by  construction)  and  proceeds 
alcnq  the  curve  toward  smaller  interatomic  distances,  using 
the  results  of  each  calculation  as  a  starting  point  for  the 
next.  This  procedure  is  particularly  important  for  charge- 
consistent  calculations  (which  is  not  the  case  here, 
however) .  The  results  cf  these  calculations  are  presented 
in  Fiqure  5  an!  summarized  in  Table  4.  The  spectroscopic 
constants  are  determined  by  a  least-squares  fit  to  the 
Fultert- Firsct f e ldcT  (1941)  modified  Morse  potential 
function  and  interpolation  is  generated  by  a  cutic  spline.  l 
The  aqreeient  with  experiment  is  quite  remarkable  -  so  much, 
in  fact,  that  one  is  a  little  reluctant  to  stake  any  claims 
as  to  its  verita tilit y. 

It  appears  that  we  have  witnessed  both  extremes  cf  cur 
earlier  expectations.  Perhaps,  thcugh,  the  latter  success 
fcr  the  potential  energy  curve  for  NO  (X2n)  can  offer 
recourse  to  a  resolution  of  the  apparent  failure  to  describe 


»Progran  written  ty  Nelson  H.   F.   Beebe,   (Quantum  Theory 
Prelect,  University  of  Florida,  Gainesville,  Florida). 


Fiqure   5.      Calculate!   and 
potential  enerqy  curves  fcr  NC(X2Tt). 


experimental   (Worse) 
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Table  <4.   Calculated  and  experimental  spectroscopic 
constants  for  nitric  cxide  (X2n)  . 


Ccnst ent 


expt 


this 
calc 


! 


calc 


r  (AO) 
e 

-1 
w  {cm   ) 
e 

-1 
w  x  (ci   ) 
e  e 

-1 
B  (cm   ) 

6 

-1 

Of  (cm   ) 
e 

D  (eV) 
e 


1.15 


1904 


1.10       1.25 


1H 


1.70 


6.5 


199  3 


20 


1.86 


€.C 


17  3  1 


2  a 


1.U5 


C.018      0.C32      0.023 


3.8 


h 


fierzterc  ( 1950)  . 
i 

Thulstrup  and  flhrn  (1972),  uses  sane  spline  fit 
as  in  present  *.crk. 


the  ionization  process.  That  is,  one  night  anticipate 
similar  success  in  confuting  the  potential  energy  curves  for 
the  varices  states  of  NO  for  a  direct  calculation  cf  the 
ionization  potentials.  Up  shall  end  the  treatment  of  nitric 
oxide  on  this  hopeful  note  since  our  real  aim  here  is 
scirewhat  different. 
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B.  Njckel  (JOC)  Surface 


Much  has  bepn  said  in  these  energy-conscious  tines  atcut 
the  iirportance  of  the  study  cf  surfaces  and  of  chemisorption 
on  them  in  the  field  of  catalysis,  and  the  recent  growth  cf 
interest  in  this  field  has  spawned  significant  advances  in 
both  theory  and  experiment.  Ultrahigh  vacuum  techniques 
coupled  kith  an  endless  variety  cf  optical,  electron  and 
other  spectroscopic  measurements  provide  the  surface 
scientist  with  the  tools  to  help  understand  the  nature  cf 
surfaces  and  the  electronic  structure  of  surf ace-adscrfcate 
complexes.1  Or  the  theoretical  side,  various  approaches  steir 
from  both  solid-state  hand  theory  and  theories  generally 
applied  tc  molecular  complexes;  however,  we  are  still  far 
from   possessing   a   reliable   computational   method   for   a 


iSee,   e.g.,    Physics  Today,   28(4), 
devoted  tc  surface  physics. 


April   1975;   issue 
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satisfactcry  description  of  the  chemisorp tion  bond.1  One 
needs  a  method  as  simple  as  extended  HUckel  theory 
(Wolfsberq  and  Helmholtz,  1952;  Anderson  and  Hoffmann,  1974) 
which  is  self-consistent  and  contains  nc  empirical 
parameters. 

The  approach  we  take  here  is,  cf  course,  a  localized  one 
with  respect  to  th*5  role  of  the  surface  in  the  chemisort picn 
bend.  In  order  to  avoid  the  effects  cf  cluster  boundaries, 
we  shall  compose  a  finite  atom-cluster  model  built  up  from 
atoms  which  retain  their  "bulk",  "surface"  and  "adsorbate" 
identities.  Fach  atcn  in  the  cluster  will  still  "feel"  the 
influence  of  at  least  its  nearest-neighbor  environment,  even 
thouqh  these  neiqhbcrs  tray  not  be  explicitly  in  the  cluster. 
This,  we  feel,  is  a  very  important  feature  *c  be 
incorporated  in  small-cluster  simulations  of  extended 
sytems.  One  would  like  to  have  a  "surface"  cluster  which 
does  not  have  surface  cf  its  own,  so  to  speak.  We  shall 
construct  our  model  in  three  steps. 

The  first  step  is  to  characterize  the  bulk,  and  by  this 
we  itean  a  self -consistent  description  of  an  atom  in  the 
metal.  In  the  case  of  nickel,  we  have  a  f ace-centered-cubic 
(FCC)  lattice  which  frca  the  viewpcint  cf  a  single  atom  is  a 


»See  Bullett  and  Cohen  (1977a)   and  references  therein  for 
a  short  synopsis  of  current  techniques. 
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nearest-neiqhtor  shell  cf  12  ether  atoms  in  a  cufccctahedral 
arranqemert  (see  top  structure  in  Figure  6),  each  of  which 
sees  the  same  environment  about  it,  and  so  on,  to  complete 
the  fully  extended  bulk  arrangement.  As  a  first 
approximation,  he  will  consider  nearest-neighbor 
interactions  cnly.  The  self-consistent  calculation  proceeds 
as  fcllovs.  We  start  with  a  free  atom  calculation1  using 
the  tyfcric  sirq  l€-dcuhle-zeta  basis  given  in  Table  5  and  a 
spherical  4-tertr  model  potential  expression  with  the 
exponents  in  Table  6.  We  then  take  13  of  these  atoms  with 
their  free-atcm  crtitals,  enerqies  and  mcdel  potentials,  and 
arranqe  +  hen>  as  in  Fiqure  6  <*of).2  The  experimental  lattice 
constant  of  4.7  a.u.  is  assumed  throughout.  The  next 
calculation  is  or  the  central  atom  in  this  cluster  with  the 
12-free-atcm  environment  taken  fully  intc  account.  Now 
comes  the  first  iirpcrtart  step  toward  self-consistency.  We 
take  the  results  of  the  central  atcm  and  place  them  at  the 
12  neighbor  sites  and  repeat  the  central  atom  calculation 
with  this  modified  e rvir enment .  This  process  is  repeated 
until  the  central  atom  and  its  neighbors  "coincide".  The 
eigenvalue  spectrum  of  this  se If -censi s*en t  "a tom-in-metal" 
is  shewn  in  Figure  7,  cenpared  to  that  of  the  free  atom. 
Cne  cculd,  in  principle,  extend  this  into  a  band  calculation 


*An  unrestricted  Grand-Cancrical   Hartree-Fock  calculation 
fcr  the  3F  state  (d8  configuration)  . 

ZFerronaqnetic  couplirq  (triplet  pairing)  of  spins. 


Fiqure    6.         Cluster   models    for 
four-fold    hole    adscrnticn    sit^. 


the  bulk,   surface  and 
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Table  5.   Easis  set  for  nickel, 


Type 


exponent 


1? 

2  "* .  3  »4 1 C 

2  s 

1 0  . 1 C  S  8 

2p 

12  .0478 

3s 

U  .  98  eo 

3p 

4.6940 

3  d 

6.7C55 

3d 

2.87  38 

Us 

2.07-51 

a? 

1.1389 

tfixed  sirqle-  ard  double-zeta  basis  cf 
Rcetti  and  dementi  (1S74). 


Table  6.   Podel  potential  exponents 
for  nickel. 


Type 


exponent 


1s 

1.5193 

1s 

3  .0366 

Is 

6.0732 

1s 

12.1464 

Even-tempered  set  (0.5,1.0,2.0,4.0) 
scaled  by  Z**  (1/3)  . 


Fiqure      7.  Calculated       phctoelect ton         spectra      of 

Ni(atotn),    Ni(netal),    and    Ni  (snr  f  ac<=>)  .        (ranqe   0-IUOeV). 
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-just  as  with  the  renormalized-atom  approach  {Hodges  et  al. , 
1972)  and  in  a.  certairly  more  elegant,  if  net  superior, 
fashion.  Such  a  venture  has  ret  been  undertaken  to  date, 
however. 


The  next  step  is  to  characterize  the  nickel  (100)  surface. 
Analogous  to  a  bulk  atom,  a  ( 100) -surf ace  atoi  has  enly  6 
nearest  neighbors  (see  riddle  structure  in  Figure  6),  if  we 
iqnore  arv  distortions  for  the  moment.  As  a  first 
approximation,  we  can  consider  the  top  layer  as  being  the 
only  layer  which  is  different  freir  these  in  the  bulk.  We 
shall  alsc,  as  in  the  bulk  case,  consider  enly  nearest- 
neiqhbor  interactions.  From  transf orm-deconvcluted  IEED 
(low  enercy  electron  difracticn)  experiments  (Landman  and 
Adams,  1974)  the  (110) -surface  layer  spacing  is  within  1$  of 
the  bulk  value  for  Ni  (FCC) r  so  no  distort iens  in  the 
qecmetry  reed  to  be  taken  into  accccnt.  Now,  since  fcur  of 
the  atoms  in  our  model  cluster  are  in  a  "bulk"  layer,  we 
held  then  fix^d  in  their  "bulk"  state  (same  orbitals, 
enerqies,  model  potentials,  etc.).  The  other  four  are 
surface  atoms  identical  to  the  central  atom.  We  will  thus 
proceed  as  before  with  a  calculation  for  the  central  atom 
takinq  irto  account  its  full  nearest-neighbor  environment 
(the  ether  fcur  surface  atoms  starting  ou*  in  their  bulk 
states) ,  and  then  iterate,  by  replacing  the  four  neighboring 
surface   atoms   with    the   results   cf   fhe    central   atou 
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calculation,  until  self-ccnsistercy  is  reached.  For 
practical  reasons  yet  to  be  worked  cut,  the  calculation  for 
the  nickel  surface  did  not  converge.  nespite  various 
damping  techniques  employed,  only  fcur  iterations  could  he 
obtained  before  strange  things  happened.  The  eigenvalue 
spectrum  for  the  fcurth  iteration  is  displayed  in  Figure  7 
alcng  with  that  of  the  free  atom  and  the  metal  atcm  results. 
It  was  also  necessary  to  held  the  ircdel  potential  fixed  to 
its  free  atcm  value,  and  for  consistency,  the  bulk 
calculation  was  redone  using  the  same  model  potential  (the 
differences  were  minor,  however).  The  shift  tc  higher 
binding  enerqies  is  scirewhat  unexpected  and  is  perhaps  only 
an  artifact  of  the  convergence  ancnaly,  but  we  shall  +ake 
what  we  can  get  and  proceed. 


The  last  step  is  the  construction  of  a  model  cluster 
representing  the  most  likely  site  for  adsorption.  Based  en 
theoretical  and  experimental  studies  cf  H,  0,  and  CC  on 
Ni  (1  CO)  surface  (Bullett,  1977b),  the  preferred  site  for 
adsorption  is  most  likely  a  four-fold  hole  site.  This  was 
also  the  site  chosen  fcr  a  previous  calculation  of  nitric 
oxide  on  nickel  by  Batra  and  Brundle  (1S76).  Nearest- 
neiqhbor  interactions  fcr  this  site  would  involve  5  nickel 
atoms  -  four  "surface"  atoms  and  one  "bulk"  atom  (see  bottom 
structure  in  Figure  6).  We  are  now  ready  tc  bring  the 
adsorbate  into  the  picture. 
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C.  Nitric  Oxide  en  Nickel  QO^)  Surface 


The  catalytic  reduction  of  nitric  oxide  (by  hydrogen) 
over  noble  metal  surfaces  is  a  reasonably  well-known  process 
(Kcbylinski  and  Taylor,  1974) .  It  is,  in  fact,  the 
underlying  process  o^  the  catalytic  converters  installed  in 
newer  emissicr-ccnt rolled  auto nobiles  (did  you  ever  wonder 
where  that  "ammonia  smell"  came  f rem?) .  We  are  not  quite 
ready  to  tackle  the  actual  chemistry  that  takes  place  en 
such  a  surface,  but  we  can  try  to  entertain  a  couple  of 
auestions  that  seme  recent  experiments  have  left  unanswered. 
Nickel  is  net  a  nctle  netal,  but  it  is  a  known  catalytic 
aqert  fcr  many  reactions  involving  the  dissociation  of  a 
nitrcqen-cxyqen  bond,  such  as  the  reduction  (by  hydrogen)  of 
nitre  cenrouncs  to  amines  (Mcrriscr.  and  Boyd,  1966)  which  is 
probably  the  most  important  synthetic  route  in  arcmatic 
chemistry  (the   introduction  of   ether  groups   into  arcmatic 
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rinqs  via  a  cia2onium  qrcup  which  is  readily  ottained  from 
primary  amines) ,  since  nitro  compounds  are  easily  prepared 
by  direct  nitration. 


So  why  dc  we  want  tc  lcck  at  nitric  cxide  on  a 
nickel(10C)  surface"?  Cne  might  think  of  this  as  a  prototype 
system  fcr  the  type  cf  reactions  just  mentioned,  but  more 
specifically,  recent  X-ray  (XPS)  and  ultraviolet  (UES) 
studies  cf  nitric  cxide  and  nitrogen  dioxide  interactions 
with  nickel  (Erucdle,  1976),  and  infrared  (IP)  data  for 
nitric  oxide  on  nickel  (Blyholder  and  Allen,  1965)  have 
posed  a  couple  cf  interesting  questions.  There  is 
reasonable  evidence  tc  suqgest  that  nitrogen  dioxide 
dissociates  or  nickel  even  at  very  lew  temperatures  (80°K) 
leavinq  ritric  cxide  and  atomic  oxygen  adsorbed  on  the 
surface.  There  is  not  enough  data,  though,  to  say  that  the 
direct  reacticn  of  nitric  cxide  with  nickel  produces  the 
sane  species  on  the  surface  (same  adsorption  site,  same 
orientation,  same  electronic  state,  etc.).  Upon  warming  (to 
3CC°K) ,  adsorbed  nitric  oxide  slowly  dissociates  into 
nitride-  and  oxide-like  species;  however,  the  former 
reaction  is  accompanied  by  the  less  o?  atomic  nitrogen. 
Subsequent  reaction  at  this  temperature  with  nitrogen 
dioxide  and  nitric  cxide,  respectively,  produce  very  similar 
XFS  and  CES  spectra,  with  less  of  nitrogen,  and,  in  the 
latter  case,   the  appearance  of  a  new1  weakly  bound  state  of 
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nitric  oxide  (easily  descrbed  en  heating  in  vacuum)  .  The 
UPS  spectra  under  these  conditions  were  un interpretable. 
There  is  qood  reascn  to  believe  that  nitric  oxide  is  irest 
strcnqly  tounc  with  nitrcqen  crierted  toward  the  surface  in 
a  four-fold  hole  site  perpendicular  to  the  surface  flare. 
However,  could  it  be  that  nitrogen  dioxide  initially 
condenses  with  cxyqen  oriented  tcward  the  surface  (since 
both  ends  have  exposed  oxygen  atoms)  and  then  dissociates  tc 
leave  oxyqen  bonded  ni  +  ric  oxide  cr.  the  surface  which  later 
dissociates  (upon  heating)  to  qive  off  nitrogen?  Could  the 
weakly  adsorbed  state  cbserved  during  extentsive  reaction 
with  nitric  oxiie  also  be  a  reversed  orientation  of  nitric 
oxide  which  then  also  dissociates  tc  yield  nitrogen?  In  the 
IR  spectrum  of  nitric  oxide  on  nickel  a  weak  unassigned  band 
occurs  which  has  been  attributed  tc  a  different  adsorption 
site.  Eut  perhaps  this  too  can  he  accounted  for  by  this 
alternate  orientation.  Thus,  if  ore  could  calculate  the 
bindinq  enerqi^s  for  nitric  oxide  on  nickel  in  bcth 
orientatiens  and  predict  the  cbserved  shifts  in  the  OPS  and 
XPS  spectra,  it  would  certainly  lend  support  to  this 
hypothesis.  A  irore  ambitious  check  would  be  to  compute  a 
potential  enerqy  curve  for  the  nitrogen-oxygen  stretch  in 
this  alternate  orientation  and  predict  the  observed 
unassiqned  frequency  in  the  IP  spectrum. 


iNew  with  respect  to  the  original  UPS  and  XPS  spectra  of 
directly  adsorbed  nitric  oxide  befcre  further  reaction  was 
allowed  tc  take  place. 
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Cur  first  calculation  was  en  a  cluster  cf  5  nickel  atoms 
arranged  as  in  Figure  7  (bottcm)  with  the  nitrogen  atoir 
oriented  tcward  the  surface  and  the  spins 
antif  errcuaqiiet  ical  ly  aliqned  with  the  ferromagnetic  surface 
(one  iright  think  of  the  entire  surface  cluster  as  a  single 
Hund's-rule  atom).  The  nickel  lattice  constant  was  left 
unaltered  at  U . 1  a.u.,  the  nitrogen-nickel  distance  was  set 
at  3.7  a.u.1  (for  all  fcur  surface  atoms),  and  the  ritrogen- 
oxyqen  cistarce  was  kept  at  its  equilibrium  molecular 
distance  cf  2.17  a.u.  As  mentioned  in  the  previous  secticnr 
the  orbitals  an  1  model  potentials  for  the  nickel  atoms  were 
fixed  at  the  tulk  ard  surface  values;  however,  the 
ccirespcr dinq  energies  were  computed  correctly  for  their 
environment.  This  is  enly  a  first  approximation  but  the 
next  level  of  approximation  would  require  (without 
reprcqr  air  ninq  a*  this  pcint)  a  much  larqer  cluster  in  order 
tc  include  a  proper  environment,  for  the  surface  rickel  atcms 
(with  prcper  stcichicnetry) .  Fcr  a  coordinate-1  site  (on 
top  cf  a  nickel  atom)  this  could  have  been  done  with  the 
nine-aton  cluster  in  Fiqure  7  (middle)  .  The  total  density 
for  nitric  oxide  (ccntr ibutio rs  from  nitrogen  and  oxygen 
atomic  crtitals)  is  shown  in  Figure  8  (bottom)  compared  tc 
molecular  nitric  oxide  (tcp) .  The  plct  plane  in  this  figure 
is   0.3   a.u.    above   the    yz-plare   to   suppress   the   1s 


lTaken  from  nickel  oxide  lattice. 


Fiqure  8.  Total  density  ccntcurs  for  ritric  oxide 
(top)  ard  nitric  cxide  on  nickel  (100)  surface  (bottom). 
(x=0.3  a  .  e.) 
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Fiqurc  9.    Siesta  density  contcurs  for  nitric  oxide  (top) 
and  nitric  oxide  or.  nickel  (100)  surface  (bottom). 
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Fiqure  10.    Fi   density  contours  for  nitric   cxide  (tec) 
and  nitric  oxide  on  nickel(100)  surface  (bottom). 


-97- 


/ 


\ 


\ 


V 

/  A 


0 


^N 


{       / 


\  (o)  ))■} 

//7Y0!\\\ 
(  I V— V  J  1 


/ 


\  \  /  / 

/  / 


98 


contributions  at  the  nuclei.  The  contours  in  these  plots 
ate  not  necessarily  at  the  same  values  so  the  comparison 
here  is  ret  too  revealinq  of  any  dissociative  character  in 
the  bond.  In  Figures  <3  and  1C,  vie  show  the  valence  sigma 
and  pi  contributions  separately  (in  the  yz-plane  this  tine) 
to  demonstrate  more  clearly  the  tendency  toward  nonbonding 
character.  The  total  energy  was  atout  1  eV  lower  than  that 
of  the  separated  nitric  oxide  molecule  and  clean  surface 
(cluster),  which  is  rather  strall  for  a  molecule  which  is 
supposedly  strongly  adsorbed  to  the  surface.  The  model 
Hamiltonian  calculation  (with  virtual  orbitals  excluded  for 
subtle  computational  reasons)  gave  valence  electron  binding 
energies  (Kccpmans*  approximation,  of  course)  which  are 
displayed  in  Fiqure  11.  The  gaussian  half-width  was  C.1  eV 
and  the  triangles  represent  a  calculation  on  the  surface 
cluster  alone.  The  inset  is  a  tracing  of  the  DPS  spectrum 
for  nitric  oxide  on  nickel  (at  8C°K) . 


The  second  calculation  was  identical  to  the  first  except 
nitrogen  and  oxygen  were  interchanged  and  the  spins  were 
coupled  in  ar  analogous  fashion  as  before.  The 
corresponding  results  are  also  shown  in  Figure  11.  Here  the 
inset  is  a  tracing  cf  the  UPS  spectra  for  adsorbed  nitrogen 
dioxide  after  extensive  reaction  at  elevated  temperature 
(3C0°K)  which  is  reported  to  be  very  similar  to  that  cf  just 
adsorbed  nitric   oxide  under   the  same   conditions.    It   is 


Figure  11.     Calculated  valence  binding  energies  for 

nitric  cxide   on  nickel  (100)    surface.    (FWHM  =  0.7   eV) . 

Trianqles   denote   surface-only   calculation.  Insets   are 
experimental  UPS  spectra  cf  Drundle  (1976). 
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under  these  conditions  where  the  alternate  orientation  is 
suspected  of  contributing  most  to  the  spectrum.  The  general 
shift  toward  hicher  binding  energies  is  evident  in  both  the 
calculated  and  experimental  results.  The  similarity  cf  the 
calculated  and  experimental  spectra  is  a  little  misleading 
unless  ere  car  nentally  substract  the  surface  contributions 
in  each  case  and  recall  that  the  experimental  "zero"  is  at 
the  Fermi  level.  F urt hermcre ,  an  attempt  to  assign  the 
nitric  oxide  peaks  was  in  disagreement  with  that  of  Batra 
and  Erundle  (197£),  as  cne  might  expect,  frcm  the  discussion 
in  the  first  section  of  this  chapter.  Ihe  XFS  spectra  is 
much  mors  definitive.  There  is  an  obvious  appearance  cf  a 
different  species  in  these  spectra  fcr  nitric  oxide  and 
ritrcqen  dioxide  under  these  same  conditions.  The  nitrogen 
and  oxygen  1s  cere  binding  energies  were  computed  using  the 
transition  state  method  where  half  an  electron  is  removed 
frcm  the  "ionized"  orbital.  The  results  are  presented  in 
Table  7.  The  ceneral  shifts  are  again  accounted  for  by 
these  two  orientations. 


As  a  final  check  on  our  supposition,  a  potential  energy 
curve  for  the  cxyger  bcrded  criertaticn  of  nitric  oxide  was 
performed.  The  oxyqen  was  kept  a^-  its  criqinal  distance 
from  the  surface  (3.7  a.u.  frcm  each  surface  atom)  while  the 
nitroqen-cxygen  distance  was  varied  perpendicular  tc  the 
surface.    The  results  are  represented   by  the  curve  ( *)   in 
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Table  7.   Calculated  and  experimental  cote  (1s)  binding 
enercies  for  ritric  oxide  en  nickel  (100)  surface. 


(eV) 


Ms 


C  1s 


Ni-N-C 

Calc 

a 
Expt 

Ni-C-N 

Calc 
t 

Expt 


H15.5 

MOC 

ui7.  a 

U03 


sua.  c 

531 

5  4  8.9 
532 


Etundle  (1976). 
t 
Attributed  tc  different  state  cf  nitric  oxide. 

Fiqure  12  compared  to  the  calculated  curve  for  nitric  oxide 
(2j)  from  the  first  section.  The  spectroscopic  constants  for 
this  curve  are  qiven  in  Table  8  along  with  the 
experimentally  observed  frequency  of  Plyholder  and  Allen 
(1S65).  The  agreement  is  mere  than  likely  fortuitous, 
censiderirq  all  the  approximations  irade  about  geometries  and 
so  forth,   but  the  implications  can  not  be  overlooked  ty  any 


The  qualitative  features  cf  cur  t wo-orientation  proposal 
are  certainly  supported  by  our  calculations.  Although  there 
is  a  number  of  ways  we   cculd  inprcve  the  calculations,   one 


Fiqure  12. 
nitric  oxide  {£) 
(*)  . 


Calculated  potential  energy  curves   for 
and  nitric  oxide  on   nickel  (100)   surface 


-1CU- 


0.2[_i  i  i  i  i  i  i  i  i  i  i  i  i  i  ii  i  i  i  i  [  i  i 


0.1 


0.0 


^  -0.1 


-0.2 


-0.3 


I  I  M  I  M  I  !  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 


NITRIC  OXIDE 


ON    NICKEL  (100) 


0 .  1 '  |  |   |   |  |  |   I   |  I   I   I  I   I   !  I   I   I   I   I   I   I   I   I   I   I   I   I   I   i   j   I   I   1   I   I  i   I   I   I   |   I   I   I   I   I  I   I   I   I 


0.0 


1.0 


2.0  3.0 

R    (BOHRI 


u.O 


5.0 


-105- 

Table    E. 

Calculat 

e<3 

spectrcsccpic   constants    for   nitric 

oxide    on    rick 

el  (100) 

sur 

face.       (oxyqen    bonded). 

Constant 

calc 

r    <AO) 

1.  19 

e 

-1 

W      (CD!         ) 

2186 

e 

-1 

W     X      (CB         ) 

1  '4 

6     € 

-1 

B     (era      ) 

1  .59 

e 

-1 

a    (cm      ) 

0.02a 

e 

C  (eV) 


-1 
E(0->  1)  (cm   ) 


7.97 


214fl  (22C5) 


Experimentally  observed  frequency  of  Blyholder  and 
Allen  (1S65)  . 


doubts  that  much  quantitative  proof  can  ever  be  presented 
for  such  a  complex  system.  We  have  not,  for  instance,  even 
considered  the  interactions  of  the  ether  itolecules  on  the 
surface  (one  in  every  hole  site  most  likely,  from  intensity 
estimates).  Temperature  and  concentration  (pressure)  were 
key  factors  in  the  experiments  and  these  effects  are 
definitely  outside  the  ranqe  of  our  current  quantum- 
theoretic  tools.  The  ultimate  qcal,  of  course,  is  to  be 
able  to   describe  the  chemistry   which  takes  place   "or"  the 
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surface  and  tc  discern  the  explicit  rcle  of  the  surface  on  a 
auantum  rcecharical  level.  As  a  first  step,  one  trust  be  able 
to  characterize  the  various  species  which  might  exist  during 
the  course  of  such  a  reaction,  for  cnly  then  can  one  begin 
to  look  at  the  dynamics  involved.  We  feel  that  a  method 
capable  cf  such  characterization  is  nearly  at  hand  and  that 
potential  enerqy  surfaces  for  such  complex  reactions  are 
-just  arccnd  the  corner.  Each  day,  more  and  ircre 
experimental  information  about  the  che mi  sorption  process 
becomes  available,  -just  as  aiore  and  more  stumbling  tlccks 
are  beinc  eliminated  in  the  quartum  mechanical  treatment. 
We  hope  that  +he  present  account  has  supplied  credible 
evidence  that  the  read  to  understanding  the  catalytic  nature 
cf  the  chemiscrption  process  on  a  quantum  mechanical  level 
is  now  being  cleared. 


CHAPTEB  X 
DISCUSSION 


The  method  we  have  tried  to  present  here  is  still  in  a 
youthful  stage  of  its  development,  having  been  conceived, 
qiven  vitality  in  the  ever-  amendable  form  of  cotrputer 
software,  and  ready  to  withstand  the  test  of  time.  Its 
maturity  will  tncst  certainly  ccme  cnly  with  growth  and 
inevitable  chanqe,  for  nearly  every  assumption  or 
approximation  that  has  b«pn  made  in  its  early  stages  can  be 
improved  upon  to  soiie  extent  if  necessary.  Before  one  can 
irake  much  constructive  (or  destructive,  for  that  matter) 
criticisu,  ttouqh,  the  rcethod  needs  considerable  testing. 
Our  primary  aim  here  was  to  present  an  application  which 
detnonstates  its  applicability  tc  a  variety  cf  systems.  Such 
a  study,  however,  dees  not  produce  any  trends  or  provide  any 
statistics  by  which  tc  judge  its  predictability  or  its 
reliability.  Studies  in  this  regard  are,  of  course, 
fcrthcomirq  and  should  prove  to  he  indicative  cf  any  need 
for  alterations  in  the  irethed. 
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It  would  te  worthwhile,  perhaps,  to  review  the 
assumptions  and  approximations  which  are  to  be  placed  under 
such  scrutiny.  The  ircst  fundairental  approximation,  outside 
the  independent-particle  model,  that  must  te  examined 
carefully  for  its  -justif ication  is  the  use  of  the  flulliken 
approximation  at  various  staqes.  It  was  th<=  use  cf  this 
approximate  form  cf  the  differential  overlap  (accurate  to 
seccnd-order  in  diatomic  overlap)  which  allows  for  the 
replacement  cf  the  constrained  eigenf unctions  by  their 
renormali2ed  unconstrained  counterparts.  And  this  is  what 
leads  tc  the  si rqle-centered  expansions.  The  next 
approximation  which  cones  to  mind  is  the  Hermitean  truncated 
(throuqh  secorid-cr^er  in  overlap)  expansion  of  the  energy- 
dependent  part  cf  the  pseudopcten tial.  The  question  here  is 
one  more  of  validity  than  of  accuracy,  insofar  as  the 
Hermiticitv  is  concerned.  The  integral  approximations  are 
perhaps  the  easiest  to  criticize,  but  these  approximations 
are  rather  straightforward  and  can  be  improved  if  the  need 
is  indicated.  The  trcd^l  Rami  ltonian  used  to  obtain 
molecular  crbitals  is  perhaps  the  nost  disputable  aspect  of 
our  treatment.  The  problem  is  not  so  much  with  the 
molecular  orbitals  or  energies  themselves  but  rather  with 
the  subseouent  population  analysis  which  on<=  would  like  tc 
employ  ir  the  self-ccnsisten t  charge  and  configuration 
procedure.  This  can  be  a  crucial  point  if  one  wishes  tc 
describe    states     which    require     heavy    mixing    of 
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ccnf iqurations.  One  might  hope  that  a  diaba+ic  single- 
conf iquration  state  in  terras  of  localized  orbitals  is  the 
closest  ircdel  of  a  real  physical  situation,  and  that  such  a 
se  If -consiste nt  procedure  would  not,  in  general,  be 
necessary.  This  subject  could  ccrtprise  a  treatise  of  its 
cwr,  so  *e  omit  any  further  discussion  here.  The  treatment 
of  spin  in  out  unrestricted  approach  can  also  be  debated  at 
qreat  lerqth  and  a  proper  discussicn  could  not  he  contained 
here.  These  last,  two  pcirts  are  quite  dependent  upon  one's 
cwn  interpretation  cf  a  real  chemical  situation  in  relation 
to  what  is  a  proper  treatment  cf  spin,  symmetry,  and  so 
forth;  and  only  further  testing  will  tear  this  out.  This 
type  of  disccurse  should  probably  be  reserved  for  more 
suitable  circumstances.  One  pcint  which  has  not  been 
mentioned  so  far  is  t he  applicability  of  the  methcd  to  large 
atoms.  This,  of  course,  entails  careful  concern  for 
relativistic  effects,  as  well  as  for  the  increased 
dimensionality.  An  attractive  recourse  here  would  be  the 
use  of  cere  pseudopotentials  which  incorporate  relativistic 
corrections  (Cas  and  wahl,  197^)  .  Such  adaptations,  as  well 
as  the  inclusion  of  various  spin  interactions  directly  into 
the  Hamiltonian,  are  currently  under  consideration. 


Tc  reiterate,  all  cf  these  pcints  demand  further  testing 
before  much  can  be  said  in  their  defense  or  otherwise,  let 
us  end   this  report  of  a   new  and  promising  approach   to  the 
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study  of  polyatomic  systems,  as  though  it  really  were  in  its 

youth  and  eaqerly   awaitinq  its  charce  to  prove   itself  in  a 

world  which   is  equally   anxious  for  the  arrival  of   such  a 
method. 
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